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Residual Entropy of Linear Polymers 
H. N. V. Temperley 


The question whether it is possible to assign a finite entropy to a polymer, bearing in 
mind the fact that it may not be a well-defined substance in thermodynamic equilibrium, is 


considered. 
that will be important at low temperatures. 


k 


It is shown that the entropy of chain configurations is the only contribution 


xisting knowledge of the thermodynamics of 


hydrocarbons enables one to arrive at a generalization of the well-known ‘flexible chain’”’ 


model of a polymer molecule. 


The new model is in good accord with experimental facts on 


the specific heat of polymers, and enables a finite entropy to be defined and to be reasonably 


accurately calculated. 


These findings are then compared with certain experimental results 


throwing light on the length of the “effective segment” of polymer chains of various types, 


this quantity being defined in the text. 


The experimental facts come from studies of poly- 


mer solutions, solid polymers, solid and liquid paraffins, and related compounds, and are 


remarkably consistent with the theoretical expectations. 


A few apparent inconsistencies and 


a few instances of alternative possible interpretations for the same facts are noted. 


1. Introduction 


In attempting to predict heats of polymerization 


from thermochemical data on long-chain hydro- | 


earbons and the monomer, it is not unusual to meet 


discrepancies between calculated and observed values | 


of the order of 1 to 10 keal/mole, Flory [5, ch. 6; 
p. 254|.' After a careful survey of the evidence, 
Flory concluded that: “The pattern presented by 
heats of polymerization of various monomers appears 
to be very well explained by two dominant factors: 
elimination of resonance energy of conjugation, and 
steric interaction between substituents.”’ In other 
words, the discrepancies can be explained purely on 
chemical considerations; the necessary readjustments 
of electronic configurations associated with the 
joining together of monomers can be much more 
mportant in some types of chain than in other types. 


However, it is sometimes pointed out that errors | 


might be caused by treating a polymer as a well- 
defined chemical substance in true thermodynamic 
equilibrium. This question may also be of interest 
in the study of other polymer reactions such as 
degradation, in which the main chain is broken up into 
fragments or greatly modified. (For reactions in 
which the main chain remains practically intact, 
errors in calculations of the configurational entropy 
will tend to balance out, since they affect both 
reactants and products.) 





The present study was | 


undertaken in an effort to estimate numerical values | 


for the residual entropy of a typical linear polymer, 
bearing in mind the four facts: first, that it is in- 
evitably a mixture of chains of many different 
lengths; second, that because of ‘freedom at the 
joints’’ any one chain may assume a wide variety of 
configurations in space; third, that in a copolymer, 
the number of ways in which two sets of monomers 
ean be lined up to form a chain can be enormously 
great; finally, that at least some of any polymer is in 
an amorphous or disordered state. It is known that 
crystallization of a polymer is never complete, and 
some polymers cannot crystallize at all. More- 
over, polymers undergo a “glass-like’’ transition at 


ED 
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temperatures at which thermal agitation is still 
considerable, so that some “frozen-in disorder” is 
to be expected. We are interpreting the glass-like 
transition as being essentially the same as that of an 
inorganic glass, certain degrees of freedom being 
put out of action as a result of a rapid rise in their 
time-constants with falling temperature. Some 
workers have claimed that there is a true phase 
transition, but a discussion of this point is hardly 
relevant in the present context. See Kauzmann 
[12], Powell and Eyring [25], Boyer and Spencer 
{2}. It might seem a hopeless task to identify the 
separate contributions to the residual entropy 
associated with each of the four factors listed above, 
but the task is very much simplified when it is 
realized that contributions of less than about R 
per mole of monomer are unlikely to be of practical 
interest. At ordinary temperatures such an entropy 
calls for a correction to the free energy of \% to 1 
keal/mole. It can very quickly be shown that the 
contributions due to 2 of the 4 causes listed above 
are negligibly small, while the entropy associated 
with copolymerization can probably be estimated 
fairly accurately in any given case. We are thus 
left only with the problem of estimating the number 
of different configurations in which a real molecule 
of polymer type can exist in a solid or liquid. It is 
this study that forms the subject of the present 
paper, but the other three possible contributions to 
the residual entropy will be considered first. 


2. Contributions From Causes Other Than 
Chain Configurations 


2.1. Communal Entropy 


We shall consider the effect of chain configura- 
tions in detail, but must first discuss a different 
effect, namely, that in a crystalline solid it is suffi- 
cient to think of each molecule as never straying 
very far from a given lattice point, but in a liquid, 
gas, or amorphous solid, it may be necessary to 
allow a typical molecule “effective access” to an 
appreciable function of the whole volume. For 








ordinary liquids this leads to the tiresome “com- 
munal entropy” problem that is still unsolved. For 
polymers, however, it is readily shown that as a 
consequence of the great size of the molecules, any 
contribution to the entropy resulting from the 
bodily transport of a whole molecule around the 
assembly is quite negligible compared with the 
contributions due to configuration effects. Physi- 
cally, the coordinates of the center of gravity of a 
chain represent just three out of the hundreds of 
coordinates that would have to be known in order 
to specify the positions of all the elements of a 
chain. More quantitatively, if we permitted the 
center of gravity of a chain to wander over the entire 
volume V of an assembly of N chains instead of 
restricting it to its own “share” V/N of that volume, 
the gain in entropy would be klnN per chain or 
Nkin N for the assembly. However, this count 
includes changes of configuration that amount to 
exchanging the positions of two or more similar 
chains, which does not lead to a physically distin- 
guishable configuration. The number of permu- 
tations of N chains among themselves is N!, so the 
proper contribution to the entropy, if the chains are 
all alike, is not more than 


Nkin N—kin (N!) = Nk 


or k per chain. It is readily shown that any inter- 
action between the chains, whether it is an attrac- 
tion or repulsion, can only result in some departure 
from complete randomness, thus making the contri- 
bution to the entropy smaller. Even if one argued 
that two polymer chains could be made physically 
distinguishable, e. g., by some slight differences in 
the substituents, so that it became physically proper 
to drop the N! term, the communal entropy would 
still only be k ln N per chain. For an assembly of 
ordinary size N~10” this would amount to some- 
thing of the order of 50 & per chain of perhaps 1,000 
monomers, whereas, as already shown above, we are 
looking for effects of the order of k per monomer 
molecule. Moreover, a preparation in which the 
number of distinguishable species of linear molecules 
approached the number of molecules present could 
only be described as an extremely bad one chemically. 
We can conclude, therefore, that the contribution due 


to communal entropy is always negligible. Besides 
effects associated with the translation of entire 


molecules, one also should consider possible contribu- 
tions due to the orientation of the molecule as a 
whole, the relative positions of all the segments 
being held fixed. The number of such possible 
orientations is certainly no greater than the number 
of physically distinguishable ways in which one 
single segment of the molecule, for example, the 
end segment, can be oriented relative to its immedi- 
ate surroundings. The order or magnitude of this 
number is readily estimated from entropy data on 
ordinary liquids and solids. 

The total entropy of an ordinary liquid or solid 
at ordinary temperatures is of the order of a few 
times /??. Even if it were assumed that all of this 
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entropy is associated with the communal effect 
together with possible orientations and “rattling 


around” of a molecule of liquid relative to its im. 
mediate surroundings, this would still mean that 


the number of physically distinguishable states of g 
typical molecule effectively excited at ordinary 
temperatures is no more than ¢* to e° (about 10 to 
100). From this it can be quite safely concluded 
that the translations and orientations of entire 
molecules of polymer of fixed configurations ean 
contribute no more than 3 to 5 times k per polymer 
chain. 


2.2. Differing Molecular Weights 


It is well known that any linear polymer prepara- 
tion contains molecules of all possible lengths from 
the monomer upwards, and it is proper to regard 
chains of two different lengths as “physically dis- 
tinguishable.’’ Since chains of more than a few 
thousand units long seldom occur in practice, the 
number of species effectively present will be small 
compared with the total number of chains. Conse- 
quently, any contribution to the entropy from this 
cause should be small. For branched as distinct 
from nearly linear polymers this argument does not 
hold. The number of paraffin isomers (C,,H,,.,) is 
known, G. Polya [24] to be of the order of 2.8", so the 
number of possible species of branched polymer with 
n=100 to 1,000 is much greater than the number of 
molecules present in an ordinary assembly, in which 
case the omission of the N! term in the communal 
entropy would be correct. Even so, as shown earlier, 
this still leads to an upper limit that is too small te 
be of significance in practice. 

The following argument from Beckett (private 
communication) enables us to set a more stringent 
upper limit in the linear case. We may suppose that 
the chain lengths are distributed about the mean 
value on the basis of the so-called “most probable” 
law of Flory [5, ch. 8, p. 319]. This may be regarded 
as the ‘‘worst possible’’ case, Jeading to the highest 
entropy of all possible distribution laws. It is 
equivalent to the assumption that chains of all 
lengths are present in proportions corresponding 
to thermodynamic equilibrium at the polymerization 
temperature. By certain methods of preparation 
the much narrower Poisson-type distribution can be 
obtained, or, by fractionation in solution, a very 
narrow spread can also be obtained. The entropy 
associated with a completely random mixture of 
chains of different species, where V, is the number 
of “‘x-mers”’ is 


S,=k In |NYIN,!! with > N,=N, 


where N is the total number of chains and M the 
total number of monomer units making up the 


assembly. Then the law of the most probable 
distribution is, Flory [5, p. 319] 
N,=N(I1- pp (2) 
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with 


N=M(1—p) 


(3 


i . oie . . . 
js the equation determining p. Substituting from 
9) and (3) into (1), approximating to the factorials 
by means of Stirling’s theorem, we find, after 
performing the summations called for by (1), 
Sy = —M{[(1—p) In (i—p)+>:p In p] (4) 
where p is related to the average length of a chain 
by (3). Thus, for a degree of polymerization 
resulting in a chain of an average length of ¢ units, 
we arrive at an entropy of ~kIn¢ per chain, or 


tint 


= per link, which again is negligible for t~1,000. 


2.3. Entropy Associated with Copolymerization 


[f two monomers A and B are present in equal 
numbers they may arranged in completely 
random fashion along each chain, unless the chemical 
kinetics of the reaction are such that A attaches 
| itself easily to B but only with difficulty to another A. 
For example, a condensation polymerization might 
be like this. In the first (random) case a residual 
entropy of k In 2 per monomer can be expected; 
in the second the chain would be almost completely 
ordered and the entropy practically zero. It should 
be noticed here that the argument must be based on 
kinetic rather than equilibrium considerations be- 
cause a chain, once formed, is not likely to break up 
again in order to assume a slightly more probable 
configuration. (This excludes conditions where 
polymer degradation is rapid.) The work is readily 
extended to the case where the two monomers are 
present in the proportions g:(l—g), the corresponding 
maximum entropy being, per link, 


be 


k(—gq In g—[l—g]In[I—g]) <k In 2 (5) 


the question whether or not it is attained again 
being outside the province of thermodynamics alone. 


3. Chain Configuration Problem 


We have now shown that with the possible excep- 
tion of the ‘copolymerization entropy,” which 
cannot be predicted solely on the basis of equilibrium 
considerations, the possible other contributions to 
the entropy are too small to be of significance. This 
study of the chain configuration problem falls into a 
humber of parts. 

(a) What can be said, a priori, about a paraffin-like 
polymer in the light of existing knowledge of the 
paraffin chain and of the random-walk problem? 

b) What can be reduced from a study of solid 
polymers? 

(c) What can be deduced from studies of polymers 
in solution? 

(d) What can be deduced from a study of related 
compounds such solid and liquid long-chain 
paraffins? 


ibs 





All of these studies have contributed something to 
the solution of the problem, and it has been possible 
to draw the following conclusions with a fair degree 
of certainty: 

(a) For a hydrocarbon-like polymer the entropy 
associated with chain configurations is no more than 
1/5 k per link. 

(b) For a rubber-like polymer the experimental 
and theoretical information is less abundant and less 
definite, but it seems safe to say that the configuration 
entropy is no more than k per monomer, 

We shall begin by asking what would be expected 
for a paraffin chain without side substituents, then 
consider briefly the possible effects of such substitu- 
tions and the presence of double bonds in the main 
chain. We shall then analyze the available experi- 
mental evidence in the light of these expectations. 


3.1. A priori Expectations 


A large number of studies have been made of the 
expected behavior of statistical assemblies of var- 
ious types of jointed chains, beginaing with the 
simplest possible case in which the “joints” are 
supposed to be entirely flexible, and in which com- 
plications due to the fact that a chain may not cross 
itself are ignored. For a survey of this work the 
reader is referred to Flory [5, ch. 10]. The idea that 
a long chain molecule can be represented in this way 
has proved extremely valuable in correlating a large 
number of experimental facts, but needs to be refined 
and generalized before it can be applied to the 
problem of calculating entropies. 

It has been pointed out by a number of writers 
that the freely jointed chain model is in many ways 
analogous to both the Langevin model of a para- 
magnetic and to the perfect gas. The analogous 
thermodynamic variables are listed in table 1. 


TABLE 1. 


Perfect gas Pressure Volume ...-..- Entropy. 


Compressibil- 


ity. 
Langevin mod-| Magnetic | Mean mag- | Susceptibility..| Entropy. 


el field netic mo- 

ment. 
Freely jointed Tension... Mean length  Young’s mod- | Entropy. 
chain of chain. ulus, 


The analogy may be carried even further in that 
all three models obey Joule’s law, the internal energy 
being entirely kinetic and independent of the exten- 
sion, volume or magnetic moment, as the case may 
be. It ean also be shown without difficulty that 
none of these models permits the use of Boltzmann’s 
principle to enable absolute calculations of entropy 
to be made. The physical reason for this is obvious; 
all three models admit of an infinite number of 
distinguishable configurations unless the chains are 
completely stretched out, the magnetic material 
magnetized to saturation, or the gas compressed to 
zero volume. Differences of entropy can, however, 
be defined, e. g., in terms of ratios of volume. Some 
of the more refiried models (e. g., van der Waals 
model of a gas and Weiss model of a magnetic ma- 
terial) which attempt to take interactions between 
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atoms into account still possess this same objection- 
able feature although they do succeed in correlating 
a large amount of experimental material. As in 
analogous problems, for which an explicit specifica- 
tion of the entropy itself is required, such as the 
discussion of vapor pressures, it becomes clear that 
quantum considerations must be introduced in order 
to get a model in which the number of distinguishable 
configurations is finite, and to obtain sensible physical 
results. It is, however, not difficult to see how such 
refinements must proceed. We are guided in the 
first place by analogy with the Brillouin treatment 
of paramagnetism, in which each atomic magnet is 
only allowed a finite number of orientations, and we 
begin by looking for some similar feature in a polymer 
chain. 

For the pure paraffin chain it is possible to state 
fairly accurately what the model must be like, 
Pitzer [20], Taylor [28]. For ethane the following 
model is consistent with both thermodynamic and 
spectroscopic evidence: the two methyl groups can 
rotate about the common single bond, but there are 
three possible energy minima separated by potential 
barriers of the order of 3 keal/mole, so that the 
coordinate @, measuring the angle of twist, is asso- 
ciated with a contribution to the energy of the form 

E=K,(1—cos 3¢). (6) 

It is natural to ask whether an energy function of 
the form (6) can be applied to each link in longer 
chain paraffins. If it can, then it becomes plain 
that the difficulty in assigning a finite entropy to 
the hydrocarbon chain disappears. For, at a suffi- 
ciently low temperature, a chain in equilibrium will 
be near one of the positions of minimum energy 
called for by (6). The effect of thermal agitation 
can then be described accurately enough as leading 
to a contribution to the potential energy of the order 


of | S \W9Or)" 


“0 Lu» ~? Where @, is the “twisting coordinate”’ 
, @ 
associated with the 7™ link in the chain. The dis- 
cussion of the problem is then that of a set of simple 
harmonic oscillators, which can be carried out 
exactly on a quantum mechanical basis. If it were 
correct to discuss a long paraffin chain in this way, 
it could be concluded that the entropy associated 
with thermal agitation should approach zero as 
T—0, leaving only the term k&ln3 per link as the 
maximum possible residual entropy, associated with 
the fact that, according to an energy-function such 
as (6), each carbon atom can be linked to the next 
in three geometrically distinguishable ways. This 
argument holds equally well even if we take account 
of degrees of freedom associated with bond lengths 
as well as bond orientations. For, at a low enough 
temperature, bond lengths must also be near their 
equilibrium values, and the effect of departures can 
also be expressed by means of a quadratic energy 
function. If we suppose further that these various 
degrees of freedom are coupled, we should obtain 
the energy that will contain product terms as well 
as squares. It is always possible, however, to find 
a linear transformation of coordinates that reduces 
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such a quadratic form to a sum of squares, at th 


same time avoiding the introduction of Produe oe 
terms into the kinetic energy. This means tha, 2 

. . necess 
the low temperature behavior of the “stretching” sco’ 
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“bending’’, and “‘twisting’’ modes of a hydrocarig, 
chain is equivalent to that of a set of simple harmonip 
oscillators, which contributes nothing to the residual 
entropy. 

On this basis it is possible to give a satisfactory 
qualitative account of the specific heats of polymer 
The degrees of freedom of independent chains shoul 
be described with reasonable accuracy by a “one. 
dimensional” spectrum of Debye type, containing 
both transverse and longitudinal modes, which would! 
lead to a specific heat proportional to T at low tem. 
peratures. The interactions between chains can bp 
described by adding to this a_ three-dimensiong| 
Debye spectrum, leading to a contribution to the 
specific heat proportional to 7°* at low temperatures 
and the phenomena of crystallization and the melting 
of crystals would also be consequences of this inter. 
action. Finally, the “glass-like transition region” 
marks the range of temperature at which the jumps 
over the potential barriers occurring in expression (6) |“! 
become frequent enough for the twisting modes of the f Kilpat 





and th 
that if 
having 
yratiol 
keal/m 
minim, 
mole. 

right | 
carbon 


chain to come into action. At lower temperatures| Pt2e? 
each chain would be “‘stiff at the joints’ but the Confir 
stretching modes would remain in operation. Ip of this 
other words, the glass-like transition ought to corre.| he%®™ 
spond to the appearance of rubber-like properties and Shepp 
to a rapid rise in the specific heat connecting low and | ‘he © 
high temperature portions of roughly one-dimensional { (U¢"“ 
Debye type, the high temperature part having a lower ability 
characteristic Debye temperature. These features plane 


seem to agree well with what is observed in practice 


he . ctee } cantly 
for the specific heats of various polymers, e. g., Furu- |“ 


kawa et al. [6], but only preliminary quantitative} i» th: 
analysis on this basis seems to have been carried out. ondlv 
It may well be that too many adjustable constants | 9¢ jt , 
are involved. The model will not be discussed further | jhe gj 
from this angle, but the above considerations suggest * perim 
strongly that it must be a reasonable first approxi- | ition | 
mation. plete, 





If an energy function such as (6) were applicable } ihe o} 
to each link of a very long chain, one could say that at ab 








the three configurations at any given joint should | youylc 
occur with very nearly equal frequency, and that this} T), 
leads to a residual entropy of the order of & per instes 
“link”. (The crude value & in 3 has to be reduced | getioy 
somewhat to allow for steric effects.) Later on We | soup 
shall assemble some experimental evidence showing } jg }jk 
that a smaller value is, in fact, more likely to be cor- | type 
rect. We are therefore prompted to ask the follow- | comp 
ing two questions: vents 
(a) To what extent is an energy function such as | 9 of 
(6) an adequate description of a long chain molecule? | orn, 
(b) How much is the entropy affected by the dele- | conc] 
tion of sterically impossible chain configurations? 1a fun 
Since a residual entropy of & per unit would still be { jhat 
noticeable thermochemically, an attempt to set & } jn o¢ 
more precise upper limit is justified. Partial answers | sider 
to both questions have been obtained by Pitzer [20] | confi 
and Taylor [28]. From a study of the thermody- | pyee’ 
namic properties of various hydrocarbons Pitzer con- | ney, 


f 





at the 
TOdue 


IS thay’ cessive bonds in the main chain lie in one plane is, 


duded that expression (6) was an oversimplification, 
gnd that the trans configuration in which any three 


ching’ 


S| nergetically, somewhat more favorable than any 
Carbon! ao 


other one, involving one or more gauche linkages. 


rMoni For brevity we shall hereafter refer to trans config- 
esidu yration of bonds as ‘‘straight’’, and to any departure 
‘from it as a “‘kink’’, although the straight configura- 
actory ion is the well-known zigzag pattern. Pitzer [20] 
es and Taylor [28] propose to replace expression (6) by 
“oul something like 

_ One. 
pt k= ky (1-cos 3¢)+&, (1-cos ¢) (7) 
V tem. 


“an be and they estimate that /, and Ey should be so chosen 
siong|\ that ua straight or trans configuration is taken as 
to the) having an energy of zero, a kinked or gauche config- 
itures | uration should have a minimum energy of about +-0.8 
velting | keal mole, and that the maxima separating the three 
“minima should have peaks of about + 4.1 keal 
mole. These figures have been shown to be nearly 
right by studies of a large number of related hydro- 
on (6 | carbons; see, for example, Pitzer [20], [21 I; Pitzer and 
of the; Kilpatrick [23]; Beckett, Pitzer, and Spitzer [1]; 
utures | Pitzer and Beckett [22]; Person and Pimentel [19]. 
t the} Confirmation that the steric energy difference is 
Ip, of this order of magnitude for butane, pentane, and 
| hexane Is provided by the spectroscopic studies of 


inter. 
710n” 
Jumps 


le 


corre- | . <- ; : 
sand | Sheppard and Szasz [26]. Such a modification in 
vy and| the energy function bas some important conse- 


: , ‘eS ») Ss , » j ‘'S ace at ’ a - 
sional | quences. It means, in the first place, that the prob 





lower | ability that any three successive bonds lie in one 
tures | Plane is not now the simple factor 4s but the signifi- 
Na cantly larger value a. zoey Where E, corresponds 
: De—F: 
ative | to the figure of 0.8 kcal/mole just mentioned. Sec- 
out. | ondly, such a chain should tend to “straighten out’’ 
ants | as it cooled, an effect that is not predicted at all by 
rther | the simpler model, but for which there is a little ex- 
sS°St | perimental evidence. Of course, the glass-like trans- 
rOxr | ition intervenes before the stra‘ghtening-out is com- 
plete, otherwise the straight configuration would be 
able } the only one that would have to be considered at all 
that | absolute zero, and the configurational entropy 
Ould | would be zero. 
this! The necessity of working with the function (7) 
per instead of (6) appears to be a consequence of inter- 
iced | actions between noncontiguous methyl and —CH,C 
We} croups, and it seems necessary to ask whether one 
ing’ is likely to have any success in applying the same 
cor | type of function to more complicated cases. Two 
OW- | complicating factors are the presence of side substit- 
'uents and the possible occurrence of double bonds 
1a or other tvpes of linkage in the main chain. The 
le? | former factor is discussed by Flory [5, ch. 6]. His 
ele- | conclusions are, in effect, equivalent to saying that 
}afunction such as (7) should still be adequate, but 
be { that the chain configuration of lowest energy is, 
+ 4) in general, no longer a simple zigzag. Steric con- 
mm | siderations imply that some other more complicated 
20) configuration may have to be chosen as the “refer- 
iy ence’” one, departures from which determine the 


energy as given by function (7) 


purposes this is not serious, but it does have to be 
borne in mind if one is attempting to compute e. g., 
the effective “spread” of a chain. In the extreme 
case of polyisobutylene only one configuration could 
be found that seemed to be sterically possible at all. 
In cases like this, the energy differences occurring 
in (7) would presumably be greater than they are 
for paraffins with a correspondingly smaller entropy 
at any given temperature. On the other hand, the 
higher potential barriers would also tend to lead to a 
higher glass-like transition temperature, so that it 
does not follow that introducing side substituents 
will automatically lower the residual entropy. 
Very little seems to be known about the theoretical 
effect of double bonds in the main chain, though it 
is generally considered that they should facilitate 
rotations of single bonds adjacent to them. 

The effect of introducing a triple bond has been 
examined by Osborne, Garner, and Yost [18] who 
studied the thermodynamic properties of dimethyl- 
acetylene CH,—C =C—CH,. They concluded that 
the potential barriers restricting relative rotation of 
the methyl groups were probably not higher than 
0.5 keal/mole, and gave reasons for thinking that 
the main effect of the triple bond is simply that the 
methyl groups are much further apart than they 
are in ethane. They argue that a double bond 
ought to have similar consequences in reducing the 
height of the barrier. The barrier in olefin hydro- 
carbons for groups adjacent to the double bond is 
about 1 to 2 keal/mole. See Kilpatrick and Pitzer 
[13], who studied ethylene, propylene, and the bu- 


tenes. This is significantly smaller than for paraffins, 
but there does not seem to be any information 
available about the existence or nonexistence of 


an effect analogous to the energy difference between 
trans and gauche configurations that were noted for 
paraffins. The indirect experimental evidence is 
certainly consistent with the idea that, in rubber- 
like polvmers, the potential barriers between differ- 
ent chain configurations are lower than in paraffins, 
and that the bias in favor of straight configurations 
may be small or nonexistent. If we could neglect 
other steric effects such as the fact that a chain may 
not cross itself our present problem would be solved. 
We could say that, as the relative probabilities of 
a given C—C bond corresponding to a straight or 


kinked configuration should be in the ratio 1:2e~#2/**, 
the partition function must contain a _ factor 
(1+ 2e-*2/*7) per link, the corresponding entropy 


per link being 


(8) 


9 
- - b 
Eo/kT _| 9) 


which falls steadily from k ln 3 to zero as E,/kT rises 
from zero to infinity. For the discussion of any 
particular type of polymer chain, it would only 
be necessary to insert the appropriate values of £, 
and of the glass-like transition temperature into (8) 
in order to arrive at the expected residual entropy. 

We must now examine the effect of removing from 


+ amie 
k | (1 tT oe 2 AEP 


For the present | the partition function terms corresponding to self- 
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crossing configurations which are implicitly included 
in the calculations leading to (8). Such problems 


are of considerable mathematical and _ physical 
interest, and although only a very few highly 


idealized ones can be solved analytically (for ex- 
ample, we do not yet know the equation of state of 
a perfect gas of rigid spheres), it seems possible to 
get an answer that is sufficiently good for present 
purposes. It is known [7, 8] that the problem of an 
actual chain cannot be treated merely as a small 
perturbation of the “unrestricted”? random walk 
problem. The effect of introducing any realistic 
type of constraint always results in a significant 
reduction in the entropy per link. Stated otherwise, 
the number of possible configurations of the restricted 
chain is only a very small fraction of that for the 
unrestricted one. Some writers have introduced 
the concept of the “equivalent random walk”’. 
They suppose that for sufficiently large N, the be- 
havior of a real chain N links long is analytically 
equivalent to that of a random walk of N/n steps, 
with a suitable choice of n and of the length of step. 
Unfortunately, recent work, Wall et al. [32, 33], 
Temperley [29] indicates that this concept is not 
strictly valid in two dimensions, nor in three dimen- 
sions for values of N small enough to be of physical 
interest, at least not if we wish to treat n as inde- 
pendent of N. Nevertheless, the concept is of value 
and will be used again later. If, for example, we 
had to do with an extremely rigid chain, strongly 
biased in favor of straight configurations, it would 
be possible to define some number, m, that measures 
the average distance between two kinks. We should 
expect this to be of the same order as the n just 
introduced, though the concepts of length of step of 
equivalent random walk and length of link of 
equivalent flexible chain are logically distinct. 

A start was made on the steric problem by Taylor 
[28] and a very similar process was used by Mon- 
troll [16] for the idealized case in which the links of 
the chain are constrained to lie along a plane square 
lattice. These authors show that the effect of re- 
moving from the partition function terms correspond- 
ing to self-crossing paths of a certain type can be 
established rigorously and very simply. The type 
of path deleted is, for the plane square lattice, any 
path containing one or more closed loops consisting 
of single squares, for example, a path passing in 
succession through the points (1,2), (2,2), (2,1), (1,1), 
(1,2). Paths containing longer loops are retained 
for the present. We summarize Montroll’s argu- 
ment for the plane square lattice, each link being 
constrained to be at right angles to those that 
precede and follow it. The first link can take any 
one of four directions. We can then proceed to the 
next link by means of either a left or a right turn, 
which we shall denote symbolically by E or F. The 
passage from the second to the third link also in- 
volves either a left or right turn, and so on. If we 
retained self-crossing paths, our generating function 
for a path n+1 links long would be simply 


(Q) 





; each term in this product, with the factors E and} 
| written in their proper order, corresponding to 4} 
possible path. We shall now show how expresgigy 
(9) is modified if we delete terms containing either of 
the successions KEE or FFF, as either of these com. 
binations corresponds to the enclosure of a singl; 
square. By induction, we can determine the gep. 
erating function S, that replaces (9) for a path 
n units long if we know it for all shorter paths. Le 
S,() and S,(F) be those parts of S, whose fing 
factors are / or F respectively. Now S,(E) may! 
end in one of the two ways FE or FEE, but in no 





others, since three /’s in a row are forbidden. Thy 
we have the relation 
S,(E)=S,_,(F)-E+S,_(F)-E-E (10 


and a similar one for S,(/). If we are concerned 
as we are here, simply with the counting of con- 
figurations in which left and right turns are equally 
weighted, & and F can be replaced in equations such 
as (9) and (10) by a single variable z, so we find 
| that S,(z) is determined by the difference equation 





S,(r) =2zS,_;(4)+2°S, (11 


o\ 0) 
which can be solved by a process quite analogous to 
that for linear differential equations, that is, by 
looking for solutions of the type S,—Ap", where A 
and p are independent of n, but may depend on z, 
Such a solution will satisfy (11) provided that 


, ' l+yi 
p—ztp—-=z 0 or 7 x ( +) (12 


giving us two independent solutions, each of which 


can be multiplied by an arbitrary function of z and | 


the two can then be added. It is readily verified 
that this is the most general solution possible, and 
the two arbitrary functions are completely deter- 
mined from the known forms of S; and S:. For the 
purpose of determining the entropy of a chain n 
steps long, as n becomes large, the answer can be 
read off directly from (12), because the complete 
solution is A,p,"+ A,p."=S, and the term involving 
p." can contribute only an insignificant amount to 
the entropy compared with that in p,". We have, 
in fact, the result that the entropy per link is reduced 
from the value & In 2 that would follow from (9) to 


— o+1] 
the distinctly smaller value & In ( : — ), as a result 


of deleting a certain type of impossible path. This 
method can readily be generalized to other types of 
lattice. For example, Taylor [28] showed that the 
removal of “pentane configurations’? reduced the 
configurational entropy of a linear paraffin from the 
value & In 3 per link that we calculate on the basis 
that successive links must be inclined at the tetra- 
hedral angle to the value k In (1++2). His result 
is equivalent to replacing equation (11) by 

2pzr—x* r(l+vy2). (12a 


pP 0 or p 
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By pentane configuration is meant any configuration 
that nearly closes upon itself in four successive steps, 
which cannot occur for any paraffin shorter than 
pentane. The fact that in such eases sizable re- 
duction in the entropy occurs when we remove a 
particular type of sterically impossible path enables 
ys to understand the failure of past attempts to 
treat the effect of steric obstruction simply as a 
small perturbation of the random walk problem. 
For a long chain, any significant reduction in the 
entropy per link means that only a very small 
fraction of the random walks gives permissible chain 
configurations. 

There is no particular difficulty in extending the 
above work to other types of lattice, including cases 
in which certain configurations of successive bonds 
are, a priori, more probable than others. Thus, in 
Montroll’s discussion of the square lattice [16] two 
adjacent links were not permitted to be in the same 
straight line, but it is quite possible to extend the 
discussion to a more general case in which right- 
angles and straights have different initial proba- 
bilities. It becomes much more difficult if we now 
attempt to remove paths involving longer loops. 
This question has been studied by the writer else- 
where [29] for the plane square lattice by methods 
differing from those outlined above, and, while the 
problem of writing down the generating function for 
all noncrossing paths is still unsolved, it is believed 
that enough is known about the properties of this 
function to enable the entropy to be determined. 


We quote the result here 
E, 2 
S=k (1 -/De-#2/*7) +. .__** 3) 
[In Ty TRT oars y2] (13) 


which, it will be observed, is of precisely the form (7) 
with 2 replaced by 2. 
the gain in energy involved in replacing a right-angle 
between two successive links by a straight. If we 
put #,=0 and compare (13) with (12) we arrive at 
the remarkable conclusion that removing just pen- 
tane configurations from random paths, the links 
inclined at the tetrahedral angle, has precisely the 
same numerical effect on the entropy as has removing 
loops of all sizes from random paths disposed along 
the square lattice, the steric effects being at least as 
severe in the paraffin chain as they are in the square 
lattice one. A removal of loops larger than those 
corresponding to pentane configurations will cer- 
tainly diminish the paraffin entropy further. We 
wish to take account of the effect of a finite /, as 
well as of steric effects, in the paraffin chain, and it 
seems safe to regard (13) as an upper limit for these 
purposes. With the value of /, of 800 cal/mole and 
a glass transition temperature of the order of 200° K, 
we should have e”2/*? so that we estimate 0.5k 
per link as an upper limit for the configurational 
entropy of a paraffin-like chain. The corresponding 
length of effective segment is six C—C links. For a 
rubber, #, may be negligible, and our estimate be- 
comes about & per link for the upper limit, and the 
corresponding effective segment is 1.7 monomer units. 


E, here has the meaning of 
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There is still one factor, not negligible, to consider, 
but there are no accurate means of estimating its 
numerical importance; this is that the occurrence of 
crystallization in polymers (in some cases up to high 
temperatures) shows that the interaction between 
chains is not negligible, and it can be foreseen that 
any steric effects between neighboring chains will 
have the effect of lowering the configurational entropy 
still further. This is probably true in a liquid 
paraffin as well as in a solid. If this effect is signifi- 
cant it would mean that the entropies actually 
occurring would be appreciably below the above 
estimates, and there is experimental evidence that 
this is sometimes the case. Various independent 
attempts to estimate the length of the effective 
segment give values for this quantity running 
appreciably higher than those that would be cal- 
culated for an isolated chain by the method outlined 
above. In this way the influence of the neglected 
factor of interaction between chains seems to be 
about comparable with that of interaction between 
distant parts of the same chain. There is some 
evidence that this neglected factor has a stronger 
influence in solids than it has in liquids or solutions, 
and this is to be expected. 


3.2. Analysis of Certain Experimental Results 


The preceding theoretical considerations suggest 
that the configurational entropy of a polymer chain 
may be significantly below the value k log 3 per link 
that one estimates from ‘elementary’ considera- 
tions. The two possible causes for this are: (a) 
The intrinsic properties of the linkages themselves, 
e. g., straight configurations may be favored by 
the nature of the interactions between neighboring 
links, as in the paraffin chain. (b) Steric effects 
tend to reduce the number of ‘‘coiled-up”’ configura- 
tions. This includes both interactions between 
neighboring chains and between distant parts of the 
same chain, while (a) is the result of short-range 
interactions within chains. What we shall try to 
estimate from experimental data is the effective 
length of segment or “average length of chain be- 
tween kinks.”’ For the paraffin chain, all the 
C—C bonds in such a segment lie in one plane and 
assume the zigzag pattern. For other types of 
chains the geometry of the configuration of lowest 
energy may be more complicated, but we shall still 
refer to it as straight, to a departure from it as a 
kink, and to the mean distance between two kinks 
as an effective segment. Strictly speaking, the last 
concept is logically not quite equivalent to that of 
length of step of equivalent random walk but the 
ratio of the two is not likely to be far from unity. 
It is, in principle, possible to distinguish experiment- 
ally whether (b) is the dominant consideration or 
whether (a) is also significant. In the latter case, 
the effective length of segment should increase as 
the temperature falls and kinks become less likely, 
and there does seem to be some evidence that this 
occurs. 

The evidence may be roughly divided as follows, 
one subsection devoted to each topic: 








1. Studies on solid polymers. Apart from the 
thermodynamic properties already referred to, as 
being in reasonable accord with the model we are 
using, there is a little further evidence giving direct 
information on the segment length in a few cases. 

2. Studies on polymer solutions. The informa- 
tion here is mainly indirect, and is usually equivalent 
to studying the effective size of the coiled up polymer 
molecule of known molecular weight. For a com- 
pletely flexible chain, the effective linear dimension 
of a coiled up molecule is proportional to N/*/, where 
N is the number of links of end-to-end length /. 
If straight configurations are so much favored that 
there is on an average only one kink per m links, 
such a chain should behave, to a first approximation, 


like a flexible one of = links. each of length Am. 


(where A is a constant, of the order of unity, depend- 
ing on the geometry of the straight configuration) 
and its mean dimension should be increased by a 
factor Am'*. We can deduce the value of m by 
comparing the dimensions of the molecule, 
inferred from experiment, with that calculated 
from bond lengths and angles on a random walk 
basis. Effects (a) and (b) both tend to increase m, 
and it is only possible to separate them in a few 


cases. 


as 


3. Studies on solid and liquid paraffins and related 
compounds. There is quite a variety of evidence 
pointing to the conclusion that such molecules, (up 
to about 30 carbon atoms), behave roughly like 
nearly rigid rods. Such a conclusion in good 
agreement with the thermodynamic properties of 
the paraffin chain discussed above. 


Is 


a. Studies on Solid Polymers 


It has been noted (and perhaps this is the most 
suggestive point) that the tensile strength of a typical 
polymer remains practically zero until a degree of 
polymerization of about 50 is reached, after which it 
begins to increase rather suddenly. It easy to 
assume that the real distinction between a “wax’’ 
and a “fiber” is that the molecules of a fiber cannot 
be sheared over one another indefinitely because the 


Is 


kinks become entangled while the molecules of 
paraffin wax are nearly straight and can slide 
smoothly past one another. One arrives at an 


effective distance between kinks of the order of 15 
to 20 units, on the basis of the rather crude assump- 
tion that molecules would have to have at least two 
kinks before they could serve as “‘double-hook”’ units, 
which could theoretically form a structure capable of 
resisting tension. 

The fact that above about a two-fold extension, 
the tension in rubber rises much more rapidly with 
the extension than the simple ‘“random-chain”’ 
theory predicts, indicates that the maximum obtain- 
able extensions, (of the order of 5 to 10 fold) really 
do represent a situation in which an appreciable 
fraction of the chains are pulled into configurations 
approaching maximum extension. This conclusion 
holds equally well whether, with Treloar [30], we 
regard the observed stress-strain curves as quali- 
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tatively explained by James and Guth’s [11] tregs. 
ment of highly extended chains, which does predie 
just such a rapid increase in elastic modulus as fylj 
extension is approached, or whether, with Flory 
[4] we attribtite this increase in elastic modulys 
mainly to the effect of crystallization. Crystalliza. 
tion is facilitated by stretching out the chains thys 
making the configurations of neighboring chains 
more nearly alike. If a 5 to 10 fold extension is 
really the maximum that can be obtained without 
rupture of an appreciable fraction of the chains, this 
would mean that 


(~).m~5 to 10(~). 
N. 


or — ~30 to 100, whereas, in an actual lightly yul- 
m A 


Mi 


canized rubber, the distances between cross-linkages 
probably correspond to chains at least several hun- 
dred units long. Thus it seems difficult to reconcile 
the observed elasticity of rubber with a value of 
1, but postulating a value of the order of 2 to § 
would improve matters considerably. Treloar [30, 
ch. 8] concluded, on the basis of studies of the double 
refraction of solid rubber under stress, that m=2, 
the precise numerical value depending on the exaet 
way in which we define m. 

A very recent study, Honnold, McCaffrey, and 
Mrowea [10] has been made using proton resonance 
methods. It is difficult to obtain actual estimates 
of segment length directly by this means, but their 
studies of silicones, and their comparison of poly- 
propylene and polypropylene oxide do indicate that 
the introduction of oxygen linkages into the main 
chain makes it more flexible. This is in agreement 
with the conclusions of Kuhn and Kuhn, |14] obtained 
from the birefringence of solutions. 


Mn 


b. 


Analysis of Data on Solutions of Polymers 


There are two basic wavs to determine the length 
of the effective segment. The most direct is to treat 
the large polymer molecule as if it were an aggregate 
of effective segments, each free move short 
distances in the solution, and assume various 
orientations. These changes in the position and 
orientation of the segments correspond to movement 
of the kinks up and down the long polymer molecules, 
the average number of kinks remaining constant. 
This movement of segments as a result of thermal 
agitation will resemble that of a solute molecule in 
an ordinary solution, except that the segment is not 
free to wander very far from the center of gravity 
of the polymer molecule. This situation is described 
bv the so-called lattice model of a polymer solution, 
as amended by Flory to cover the case of extremely 
dilute solutions. See Flory [5, ch. 12] for a summary. 

Using this treatment, we can estimate the size of 
the effective segment in several ways. 

(a) From the osmotic pressure. Since this is an 
equilibrium property, we can use solution theory 
directly. The limitation is set by the availability 
of data. The solutions have to be very dilute for 
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theory to be acceptable, whereupon the 


existing 
pressures to be measured become rather low for 
accurate measurement. 

(b) From the viscosity. 
motion of a segment through the solvent is associated 
with the displacement of an equivalent volume of 


The idea here is that any 


solvent molecules. These molecules have to be 
removed from the liquid, leaving a “hole” big enough 
to contain the segment. The segment then enters 
this hole and the space it leaves vacant in the liquid 
is filled by the displaced molecules. Independently of 
any particular model of the liquid, we can argue that 
the energy barrier associated with such a process 
must be of the order of magnitude /pr, where / is the 
latent heat of evaporation, p is the density of the 
solvent and ¢ the volume of the effective segment. 
Thus the probability of such a process may be 
expected to be proportional to exp (—/pr/kT), and 
the viscosity of the solution should contain a factor 


that is the inverse of this. If so, the activation 
energy /pr can be estimated by means of the 
customary plot of log (viscosity) against 1/7. For 


a survey of the conclusions from these two methods 
we quote Powell and Eyring [25, p. 191]. “It is 


worth noting that the segment lengths estimated 
from Osmotic pressure data (18 to 33 atoms for 
polystyrene and polyethylene oxide) and from 


melting point data (20 to 25 atoms for hydrocarbons) 
are in agreement with the segment lengths from 
viscosity data (20 to 25 atoms for hydrocarbons, 
28 to 34 atoms for polvesters).”” The reader 
referred to this review article for a survey of 
work leading to these conclusions. 
Eyring and his collaborators have attempted an 
alternative estimate of segment length from viscosity 
data, on a viscosity formula of the type 


is 


the 


based 
n=A exp (u/kT), 


in which attempts are made to estimate A, as 
well as uw, from the experimental results, and to 
compute a theoretical value of A in terms of segment 
length. In practice, this method fails to give useful 
results, but this failure does not mean that the esti- 
mates of uw are invalid also. Suppose, for example, 
that w is varving slightly with temperature, say 


according to the law u=w(1+a T). The customary 
logarithmic plot of the viscosity data against 1/7 
would now vield uw and Ae*“e* in place of uw and 


and a value of a of no more than a few tenths of one 
percent per degree could lead to a completely false 
experimental estimate of A. A similar difficulty 
occurs in the analysis of thermionic data in metals. 

(c) A third method of estimating segment lengths, 
due to Kuhn and Kuhn [14], rests on the study of the 
double refraction of a polymer solution that is being | 
steadily sheared, the shearing flow tending to orient 
the typical segment and stretch out the polymer 
molecule in a manner that depends essentially on the 
length of the typical segment. By a rather compli- 
cated method of analysis, these authors arrive at 
results ranging from segment lengths of the order of 
30 carbon-carbon links for paraffin-like polymers 
ranging down to 1 to 3 monomer units for rubber, 
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hydrogenated rubber, and polymers which contain 

-O linkages in the main chain. Such linkages, by 
increasing the spatial separation of the —CH,C— 
groups, no doubt make the chain more flexible, as 
already suggested above. 

The second basic method of attack on this question 

based on an attempt to determine the linear 
dimensions of the molecules in solution, and to 
examine how these differ from the values to be 
expected on the completely random-walk basis. This 
should give a measure of m, since we have pointed 
out that the ratio, actual dimension to calculated 
dimension, should be approximately equal to m?. 
The actual dimension is obtained by an ingenious 
analysis of data on the viscosity of dilute solutions, 
due to Flory [5, ch. 14]. This new method of analysis 
of viscosity bears no relation to the method we have 
already described, but is not inconsistent with it. 

The very dilute polymer solution can be regarded 
as a mixture of roughly spherical aggregates repre- 
senting the molec ules of polymer, separated by fairly 
large regions of pure solvent. Flory’s idea is to 
suppose each molecule of polymer replaced by an 
equivalent solid sphere, and to attribute the excess 
of the viscosity of the solution over that of the pure 
solvent to the fact that the solvent cannot flow freely 
through the polymer molecules, but is partly diverted 
around them. One can then apply Einstein’s well- 
known treatment of the viscosity of a liquid contain- 
ing a random distribution of solid spheres. Thus, 
according to this approach, we consider that a whole 
polymer molecule obstructs the flow of the solvent 
to an extent that depends on its spread in space, 
and the motion of individual segments not 
considered. 

Flory shows that, if the temperature is properly 
chosen in relation to the solvent and concentration, 
the spatial distribution of polymer molecules in the 
solution is indeed random; and further, that the ex- 
tension of the molecule in space is then very nearly the 
same as what we should calculate on the assumption 
that the interactions of solvent molecules and 
polymer segments with one another and among them- 
selves are ‘equivale nt. At this temperature, called 
by Flory the © point, the solvent becomes, so to 
speak, ‘‘neutral’’, any actual differences in these 
interaction energies being approximately balanced 
by the effect of thermal agitation. Various other 
definitions of the @ point are possible; for e xample, it 
can be shown that it should correspond to the limit 
of solubility of a sample of polymer of extremely 
high molecular weight. A solution at its @ point 
should then contain a fake ally random spatial 
distribution of molecules in the solvent, and the 
average chain configuration of these molecules should 
be that which we should predict on the basis of 
neglect of long-range interactions between segments 
(apart from steric effec ts). 

Flory’s analysis [5, ch. 14, p. 618] leads him to con- 
clude that “the observed end-to-end distances, 
unperturbed by long-range interactions, are con- 
sistently greater than would hold if bond rotation 
were complete ‘ly unrestricted. Steric effects of one 
sort or another probably are responsible.” The 


is 
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factor measuring the “expansion of linear dimen- 
sions’’, which should be a measure of m'?, ranges 
from 1.60 for gutta-percha to 4.9 for cellulose tribu- 
tyrate, in satisfactory qualitative agreement with 
figures for m already calculated for paraffin-like and 
rubber-like polymers by other means. 

It is of some interest to ask what is likely to be the 
relative importance in causing this expansion of the 
molecule of the two distinct factors met in discussing 
the hydrocarbon chain, the fact that two carbon 
atoms may not be in the same place, and the fact 
that there may be correlation between the directions 
of e. g., neighboring (—C bonds. The first of these 
factors should lead to an expansion that remains 
the same at all temperatures, while the second 
should lead to an expansion that increases as the 
temperature falls. Calculations on a_paraffin-like 
chain, Taylor [28], Kuhn [15] indicate that the end- 
to-end length should increase by a factor of 1.2 as 
the temperature falls from + 100° C to 0° C. Flory 
[5] finds evidence of just such an effect for various 
solutions of cellulose tributyrate, each one studied 
at its 0 point. This is also the very molecule that 
shows the biggest departure from the random-walk 
configuration. 

On the basis of all this evidence, we seem justified 
in concluding that: 

a. The absence of configurations in which two 
carbon atoms coincide increases the linear dimen- 
sions by a factor of the order of 1 to 5 above what 
they would be on a random-walk treatment. This 
effect seems to be the only important one in rubber- 
like materials. It implies that only about half the 
linkages between isoprene units can be regarded as 
uncorrelated with one another, and leads to a residual 
entropy of '; & In 3 per isoprene unit. 

b. For a paraffin-like chain there is often an addi- 
tional effect, due to correlations between neighboring 
links probably attributable to interactions between 
the CH,—C groups in the main chain. This seems 
to be capable of raising the effective value of m from 
a value of the order of 2 that seems to follow on the 
basis of factor (a) alone to a value of the order of 20 
or more. The corresponding residual entropy is 
then of the order of 0.1 to 0.2 k per monomer mole- 
cule. It would not be correct to quote a figure as 
low as 0.05 k corresponding to an effective segment 
20 units long, because of the presence of the second 
term in eq (8). 

These values agree remarkably well with thos 
based on the purely theoretical study, bearing in 
mind that one factor, the interaction of the chains 
in a solid with one another, was not considered. 
It is satisfactory that the apparent values of the 
length of the effective segment are, for both rubber 
and the paraffin-like polymers, a factor of 2 or so 
greater than we estimate theoretically for an iso- 
lated chain. These experimental studies are almost 
all based on dubious assumptions, or depend on 
numerical constants that have been evaluated very 
tentatively. Any one of them, taken alone, would 
not be very impressive, as comparatively small 
errors in some of the numerical constants might 


make a considerable difference to the calculated 
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values of m. The results of a number of quite 
independent investigations all give calculated values 
of the same order of magnitude, and their cumula- 
tive value as evidence is very considerable. 

We shall close with a survey of the properties of 
liquid and solid paraffins and similar bodies of chain 
lengths between about 10 and 60, which seems to 
confirm the above conclusions obtained directly 
from polymer studies. " 


c. Analysis of Data on Paraffins and Allied Series 


Convincing evidence that paraffin molecules up to 
about 40 units long pack like rods in the solid state 
is provided by the X-ray work of A. Miiller. For 
an account of his work the reader is referred to a 
review article by Daniel [3]. In addition to deter- 
mining the structure of the ordinary solids, Miiller 
was able to elucidate the nature of a transition that 
occurs a little below the melting point for paraffins 
of lengths between 16 and 32 C atoms. He was 
able to show that this transition is associated with 
a change in crystal structure, the packing changing 
from a type natural for flat ribbons to one appro- 
priate to practically symmetric round rods. 

The following evidence that in the high tempera- 
ture state the molecules really are rotating as rigid 
rods rather than twisting is available. In the first 
place, the entropy change associated with the transi- 
tion has been found to be practically independent of 
chain length, which is what we should expect if the 
entropy change is mainly associated with a change 
in the two-dimensional packing of rigid rods. But 
if only short segments of the chain were rotating at 
any one time and some of the bonds were being 
twisted, one would expect the entropy of transition 
to be proportional to the number of segments, and 
thus to the length of the chain. Also, the transition 
temperature is a linear function of the chain length 
y, of the form A+By. This, combined with the 
result that the entropy of the transition is practically 
independent of chain length, implies that the energy 
change associated with the transition must be a 
similar linear function. This is just what would be 
expected if the chains are behaving as rigid rods, 
one contribution to the energy coming from the 
interaction between any chain and those lying along- 
side it, which should be proportional to the length 
of the chain, y. A second contribution might be 
expected from the interactions of methyl groups 
belonging to end-to-end neighbors, which should be 
practically independent of chain length. 

This type of transition occurs only over a limited 
range of chain lengths simply because, at the top 
and bottom of this range, the paraffin melts before 
undergoing the transition. The melting tempera- 
tures of various homologous series are well repre- 
sented by the empirical relation 

Tm = (C+ Dy) /(E+ Fy) (14) 
in contrast with the linear relation for the transition 
temperature just mentioned. It is now natural to 
ask whether eq (14) is also to be regarded as evidence 
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hat paraffin chains are behaving as practically rigid 
ods. In view of the tentative nature of present 
theories of liquid structure, any attempt to be dog- 
matic on this point would be out of place. It does, 
however, Seem worth while to mention two possible 
interpretations of (14), which, although incompatible 
with one another, are both in agreement with the 
picture of a fairly long effective segment. It could 
be argued that since for all homologous aliphatic 
series, a relation such as (14) seems to hold, and that 
the melting temperature settles down to a practically 
constant value for m=30, that this represents the 
length of the effective segment; that lower hydrocar- 
hons are behaving like rigid rods in both solid and 
liquid states, while higher ones are behaving as if 
they were composed of effective segments of about 
0 to 30 links long, each of which is capable of moving 
ground in the liquid as a unit, corresponding to the 
motion of individual kinks up and down the chains. 
We have already met this picture in discussing 
polymer solutions. On the other hand, an interpre- 
tation of relation (14) can be given on the basis that 
the chains behave like nearly rigid rods for all values 
of y for which reasonably pure preparations can be 
made. The energy change associated with melting 
should be the sum of interactions between side-to- 
side and end-to-end neighboring chains, and_ this 
should be a linear function of y, as before. In the 
liquid near its freezing point it may be supposed 
that the chains are still packing in a roughly two- 
dimensional fashion, but somewhat more loosely 
than in the solid, and that each chain also acquires 
some freedom to move lengthways as well as side- 
ways. Such a picture leads to an entropy change 
consisting of a constant term, associated with the 
possibility of sideways movement and rotation, that 
we have already met in the solid-solid transition, 
with the addition of a term, roughly proportional to 
chain length, associated with a possible lengthways 
sliding of the chains, which would be absent in the 
solid transition. Such behavior of the entropy and 
energy of transition leads directly to a relation like 
14). For a review of the experimental material on 
which these considerations are based, and also for 
some extremely interesting material on esters, which 
s also in reasonable accord with the “rigid-rod”’ 
picture, the reader is referred to Daniel’s review 
article [3]. Simanouti and Mizushima, |27] from a 
spectroscopic study of liquid cetane (C,;H3.), con- 
cluded that this liquid was probably a mixture of 
the various possible forms of molecule containing 
just one kink or gauche configuration, all the re- 
} maining linkages being in the trans configuration but 
the molecules in the solid being entirely of the 
all trans (coplanar) type. Thus, the effect of the 
interactions between neighboring molecules seems to 
be somewhat stronger in the solid, a result that would 
be expected and this is also in keeping with the idea 
of an effective segment of ~20 units. 
We conclude by a brief reference to the work of 
Moore, Gibbs, and Eyring [17], who succeed in inter- 
preting the temperature dependence of the viscosity 
of liquid paraffins on this same basis of approxi- 
mately two-dimensional packing, so that the flow of 











such a liquid resembles the distortion of a crystal 
under stress by the sliding of layers of molecules over 
one another, rather than the independent motion of 
spherical molecules characteristic of ordinary liquids. 
The latter is what we should expect of the paraffins 
if there were more than 2 or 3 kinks per molecule, 
because then they would approximate to spheres far 
better than to rods. 


4. Discussion and Conclusions 


This enquiry has led into a variety of different 
fields and provided us with evidence on the chain 
configuration problem of very varying degrees of 
relevance and probable reliability, but the over-all 
consistency of the experimental evidence is surpris- 
ingly good, as is the consistency with the theoretical 
expectations. We have already met a number of 
cases in which the same data might be interpreted in 
more than one way, and, in the present state of our 
knowledge, both ways may well be wrong. Some of 
the investigations involve numerical constants evalu- 
ated by dubious methods, but, for the purposes of 
the present enquiry, the fact that these investiga- 
tions were carried out quite independently of one 
another is a powerful argument in favor of the quali- 
tative conclusions. It is right to mention at this 
point a few apparent inconsistencies which cannot be 
explained away so easily. The general picture is 
that a paraffin-like chain is fairly rigid, and that the 
introduction of substituents only tends to increase 
the rigidity, hence a residual entropy of about 0.2 k 
per link, but that the introduction of a double bond 
or a C—O linkage into the main chain should intro- 
duce flexibility, and thus put up the residual entropy 
to as much as 0.5 to 1.0 k& per link. And yet, (a) ac- 
cording to theory, and also according to Kuhn and 
Kuhn [14] and Honnold, McCaffrey, and Mrowca 
|10] the polyethylene oxide chain should be very 
flexible, whereas the behavior of this solid is not 
rubber-like, and osmotic pressure data, Powell and 
Evring [25], indicate an effective segment of the order 
of 20 to 30 atoms. (b) Kuhn and Kuhn [14] studied 
a hydrogenated rubber, in which all the double bonds 
were saturated. This should behave like a hydro- 
carbon chain with one substituent methyl group per 
isoprene unit, and should therefore form a stiff chain, 
whereas, in solution, it appears to be at least as 
flexible as ordinary rubber. 


5. Possible Further Work 


The main conclusions seem confirmed by a number 
of independent pieces of evidence, but the investiga- 
tion does suggest possible further lines of enquiry. 

(a) The idea of a possible numerical relationship 
between effective length of segment and fraction of 
material crystallized seems likely to throw light on 
the crystallization of rubbers, if it could be made 
more precise. 

(b) It might be worth studying nonsaturated hy- 
drocarbon liquids by X-ray and viscosity methods in 
order to see whether such molecules pack as approxi- 
mate spheres, as they should if they contain a number 
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of fairly flexible bonds, or whether they pack as long 
rods. One liquid that might be interesting to study 
in these ways is squalene, since it consists of six 
isoprene units, and the average molecule should 
contain 2 or 3 kinks. 
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Molecular Weights of Thermally Degraded 


Polymethyl Methacrylate 


Victor E. Hart 


Molecular weights of polymer residues from the thermal degradation in vacuum of two 
polymethyl methacrylate samples, one polymerized with benzoyl peroxide and the other 


iurnal of Research of the National Bureau of Standards 
without anv initiator, 
constant temperatures. 
of rate of volatilization and rate of change of 1/M,, 
1/M, R 
molecular weight 
mechanism in which intermolecular chain 
the primary, 
1. Introduction 
Polymethyl methacrylate is of particular interest 


in the study of polymer degradation because the 
products volatilized as a result of the thermal de- 
composition in vacuum consist almost entirely of 
the monomer, and because the molecular 
weights of the residues can be conveniently measured. 

Relatively little has been reported in the literature 
| to 6]! with regard to the degradation of poly- 
methyl methacrylate prepared without the use of 
any initiator. Madorsky [1] and Straus’ and 
Madorsky [2] found such a polymethyl! methacrylate 
to be considerably more stable (in terms of volatili- 
zation rates) at elevated temperatures in vacuum 
than a sample polymerized with the aid of benzoyl 
peroxide. They found that in the course of 
pyrolysis, the residue the benzoyl peroxide- 
polymerized polymer eventually achieved a stability 
after about 50% of volatilization) comparable to 
that of the polymer polymerized without any 
initiator 

The main object of the present work was to deter- 
mine molecular weights of the above two polymethy! 
methacrylate samples as functions of extent of vola- 
tilization at constant temperature for thermal deg- 
radation in vacuum. Also, it was desired to elucidate 
the nature of the temperature effect. In the case 
of the benzoyl peroxide-polvmerized sample, it was 
of particular interest to obtain data beyond the 
first 50 percent of volatilization, where complicating 
effects of benzoyl peroxide might no longer be present. 


also 


also 


of 


2. Polymethyl Methacrylate Samples 


A was analvzed and found to contain 
approximately 0.6 percent of benzoyl peroxide. Its 
weight-average molecular weight was 146,000.2. The 
monomer used in the preparation of sample B was 
degassed and sealed in an evacuated tube without 
any initiator. It was allowed to polymerize in bulk 
at —2 C. The weight-average molecular 
weight of the resulting polymer was 5,030,000 (see 
footnote 2 


Sample 


5° to —35 


! Figures in brackets indicate t 
The molecular weights 
using the light-scattering method 


terature relerences 
re determined by Louis 


Williams of the Bureau 


it the end of this paper. 


have been obtained over wide extents of degradation at 


ratio of residue weight to original specimen weight, and M,, 
The results are explained on the basis of a free-radical chain-reaction 
transfer ( 
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various 


A linear relation, dependent upon the temperature but independent 


is found when log RF is plotted against 
weight-average 
is 


activation energy=29 kilocalories) 


or perhaps the sole, cause of the molecular-weight decrease 


3. Experimental Methods 


The vacuum apparatus used in pyrolysis was that 
described previously [7] for the study of the rates of 
volatilization of polytetrafluoroethylene by the gas- 
pressure method. The specimen was inserted into 
a weighed Pyrex glass tube, 45 mm long and 6 mm 
inside diameter with one end closed, and was weighed 
on a semimicrobalance. The pyrolysis apparatus, 
with specimen in place, was kept under vacuum over- 
night to minimize desorption of gases from the walls 
of the apparatus during pyrolysis. The cylindrical 
external heater was so constructed that an essentially 
uniform temperature could be maintained over the 
entire specimen. The heater, after being kept at the 
desired temperature overnight, was moved into po- 
sition in two steps, 10 min apart, during which time 
the current was increased by a_ predetermined 
amount and duration. The specimen reached the 
desired temperature in about 2 min after the heater 
was placed in final position. The temperature was 
kept within +0.5° C. The volatile products were 
condensed during pyrolysis in a liquid-nitrogen trap 
located between the specimen and the mercury- 
vapor pump. ' 

The polymer residue was dissolved out of the Py- 
rex glass tube with a minimum amount of distilled 
ACS-grade chloroform by using the small reflux 
apparatus shown in figure 1. Moderate reflux rates 
were used to avoid foaming, and complete solution 
usually occurred in less than 2 hrs. The Pyrex tube 
was then weighed to make certain that the residue 
was completely dissolved. The solution was al- 
lowed to stand overnight prior to viscosity measure- 
ment. 

One milliliter of solution was required for an in- 
trinsic viscosity determination. he semimicro- 
dilution viscometer and the techniques used have 
been described by the author in a previous publica- 
tion [8]. Kinetic-energy corrections and errors due 
to drainage, surface tension, solvent evaporation, 
and the effects of shear rate of the polymethyl 
methacrylate solutions were negligible with this vis- 
cometer. The temperature of the water bath was 
maintained at 27.30 +0.01°C, The flow time of the 
solvent was 195.7 sec. The latter was redetermined 
after each intrinsic-viscosity measurement, and at 
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FiGuRE 1, Dissolving fused polymer from a lube 
A, Weighing bottle containing solvent; B, tube containing 
fused polymer; C, indentation in the glass; D, tapered ground 
glass connection; E, condenser tube; F, beaker for keeping out 


dust; G, heating bath 


no time was it found to have changed. Average 
deviation of flow time for solution reruns was less 
than 0.1 percent, and at least one rerun was always 
made. 

A graphic method described previously by the 
author [9] was used for accurately calculating 
intrinsic viscosity from a single specilic-viscosity de- 
termination. This method is applicable to a wide 
range of specific viscosities. However, for conven- 
lence in measuring the speciic viscosity, a value close 
to 1.0 was generally obtained by adjusting the size 
of the original polymer specimen and the concentra- 
tion of the solution according to predictions based 
upon the trends of the data. 


4. Results 


The weight-average molecular weight, 1/,, was 
calculated from the intrinsic viscosity, !n], using the 
equation [n]=3.54*107°M,°"*". The constants in 
this equation were calculated from intrinsic-viscosity 
and light-scattering molecular-weight measurements 
of the original polymethyl methacrylate samples A 
and B. ‘the above equation agrees satisfactorily 
with those reported in the literature for polymethyl 
methacrylate in chloroform [10, 11]. 

In table 1 are listed molecular weights of residues 
from pyrolysis of the two polymethyl methacrylate 
samples. Prior to pyrolysis at 320°, each specimen 
of sample A was preheated by raising it in a con- 
sistent manner from room temperature to 260° in 
55 min in the vacuum apparatus. This prevented 
excessively rapid initial volatilization at 320° C 

In figure 2, weight-average molecular weight is 
plotted in terms of percentage of the original molec- 
ular weight against percentage volatilized at 320 
C and at 150° C. Because of a significant tempera- 
ture effect, it would be incorrect to join the two 
curves for sample A. 
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The nature of the molecular-weight change is 


better illustrated in figure 3, where log R is plotted 
against 1/.\/,, R being the ratio of the residue weight 
to the weight of the original specimen. Straight 
lines result, the slopes of which depend upon the 
temperature. (The first 48% of the volatilization 
of sample A in the experiments at 320° C occurred 
at a high rate, before the temperature of 320° C could 
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be reached. For this reason it was experimentally 
impracticable to determine the molecular weights 
within this range of volatilization at this tempera- 
ture.) From figure 3 it may be inferred that in the 
type of plot shown in figure 2, the higher the mole- 
pat weight, the more abruptly the molecular- 
weight curve drops, and also the higher the tempera- 
ture, the more rapid the drop. 


5. Theory for Log Ff as a Function of 1/M, 


m/m>, Where m= weight of specimen 


Defining: R 


remaining at time, ¢, and m,=original weight of 
specimen; A,=(—100/R)(dR/dt) 230.3 (d log 
R\/dt=rate of volatilization in percentage of the 


residue per minute; A,—d(1/M,,)/dt=rate of change 


of reciprocal weight-average molecular weight. As a 
consequence of the definitions, 
> 
K: 930.3 d log R (1) 
4 da M,.) 


The object of the ensuing treatment is to indicate 
why K,/K; might be expected to be very nearly 
constant for a given py roly sis temperature. 

A free-radical chain-reaction mechanism (sta- 
tionary type) is assumed, in which the entire drop 
in molecular weight during pyrolysis is due to an 
intermolecular chain-transfer process. The initia- 
tion step of the chain reaction may be ignored, as 
no kinetic expression for this step is required for the 
present purpose. 

The volatilization products of the thermal decom- 
position of polymethyl methacrylate consist almost 
entirely of the monomer [2]. As is customary in 
this case, it is assumed that each polymeric free 
radical resulting from the initiation process decreases 
in size by the consecutive splitting off of monomer 
units from its free-radical end (depropagation). 
Svymbolically, for kinetic analysis 


h, 
P,P, +M1 
ke 
P, —P, ,+Mt 
ete 


where P; is a polymeric free radical of degree of 
polymerization, 6, and \/ represents a monomer 
molecule. 

If A, is the number of monomer units splitting off 
per gram of residue per minute, 


Kyu=kd Ps}, (2) 
where |/’-] is the number of free radicals per gram of 
residue, and ky is the depropagation-rate constant. 
Since this depropagation step contributes almost the 
entire amount of the volatilization products, the rate 
of volatilization, Ay, may be expressed similarly: 


100.1100 
N 


k, ; (3) 


talP-] 


where NV is the Avogadro number. (The molecular 
weight of the monomer is 100.1.) 

The chain-transfer mechanism is assumed to be of 
the usual hydrogen-transfer type, in which a free- 
radical end acquires a hydrogen atom from a carbon 
atom in the backbone of another chain and becomes 
a polymer molecule with a saturated end, S. The 
molecule attacked splits into two fragments, one of 
which forms an unsaturated end, UU’, and the other 
fragment contains a free-radical end. The latter 
polymer fragment commences to split off monomer. 
The transfer process may be represented for kinetic 
analysis as follows (intermolecular type): 


ir 


I 
P,_,+P,HP,—P,_,S+P;P, 


b-2z 


P: P,->P;+P,U (rapid step), 


where H is the transfer hydrogen, z is the number of 
monomer units split off from the free radical prior 
to the transfer act, and P, and P, represent the two 
polymeric fragments formed from the molecule 
attacked at random, P,HP,. 

If K,, is the number of transfer acts occurring per 
gram of residue per minute, &,, is the transfer-rate 
constant, and [H] is the number of points at which 
transfer can occur per gram of residue, which is 
constant, then 

K,,=k,|P-)(H). (4) 

Because intermolecular chain transfer is assumed 
to be the sole cause of the molecular-weight decrease, 
it follows, as a close approximation within the ex- 
perimental ranges of volatilization at 300° to 340° 
C, that 


df 


N 
i.) 


dt (5) 


K ws, 


where MM, is the number-average molecular weight. 
The following is also assumed as a good approxima- 
tion within these experimental ranges [12]: 


A 
—“—-q=constant. 6) 
uM, ! ( 
Substituting for VM, and K,, in eq (5), 
, q > ny 
K, +f helt -}{H]}. (7) 
N 
Termination of a reaction chain is assumed to 


occur merely when a free radical becomes reduced 
by depropagation to monomer size or to the point 
where it escapes as a polymeric fragment along with 
the monomer being evolved. This may be called 
“termination by free-radical depletion.” 
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mains _ 8) 
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It has been assumed that g and [H] are constant. 
If it is further assumed that the rate constants, ky 
and k,,, are indeed constant at a fixed temperature, 
it follows that K,/K, should be constant, and refer- 
ring to eq (1), that (d log R) [d(i M,,)| should be 
constant. It is clear that the theory predicts for a 
given sample that a plot of log R against 1/\/, 
(as in fig. 3) will be a straight line for pyrolysis at 
a fixed temperature, independent of change in mag- 
nitude of A, and A;, and that the slope of the plot 
will vary with the temperature, depending on the 
ratio, ke/k;y-. 

For the purpose of testing the freedom of the experi- 
mental slopes from dependence upon the magnitude 
of the rates, A, and K;, it is fortunate that these rates 
have been found to vary with extent of degradation 
at 320° C. Both samples A and B showed decreasing 
rates, and at the end of the experimental ranges of 
pyrolysis the rates were about a third of the initial 
values in both cases. These results further support 
intermolecular chain transfer as the primary, or 
perhaps the sole, cause of the molecular-weight 
change. The experimental rate curves covering wide 
extents of degradation at 320° C, together with 
studies of the effects of benzoyl peroxide and of in- 
hibitors, are planned for inclusion in a separate pub- 
lication, in which the initiation mechanism may be 
elucidated. 

In the case of sample A, 48 percent of volatilization 
occurs before the residue can reach 320 If the 
actual values of log R versus 1/M,, for those specimens 
of A for which data were obtained at 320° were 
plotted for the first 48 percent of volatilization, the 
curve would reflect the lower temperatures in this 





range and the corresponding changes in_ slope, {; 
Therefore, the straight line obtained at 320° does not |} 
extrapolate back to the initial molecular weight gt 
log R=0. Sample B incurs about 8 percent of vola- 
tilization before it can reach 320°. However, this 
amount is quite small for initial curvature to be 
indicated, and the straight lines do appear in this 
case to extrapolate to the initial molecular weight } 
at log R=0. 

The slope for sample B is 1.67 times that for sample 
A at 320° C. 
able at present for this difference. However, some 
heterogeneity in the molecular structure is  sug- 
gested in the case of sample A, at least after the 
temperature of the sample has been raised sufficiently 
high. A possible reaction of some sort involving 
benzoyl peroxide and groups along the polymer | 
chain, and resulting in a few points at which the 
probability of transfer is greatly enhanced, may be 
considered. , 





No satisfactory explanation is avail. | 


6. Activation Energies 


Rearranging eq (8), 


1.001 
g{H) 


10% Kk, 


kee = 
, K 
The constants may be combined and replaced with 
Cs: 

Kk, 


kp=C, kK 


f 


(9) 


Taking the logarithm of the latter expression and 


then the derivative with respect to 1/7, where T is 
the absolute temperature, 
dlogk, dlog K, dlog Kk, , dlogk ' 
d/T)~ dQ/T)~ da/T) *da/Ty “ 
It follows that 
E,~=E;,— E+ K, 1] 


where /,, and ky are the activation energies for the 
chain transfer act and the depropagation act, re- 
spectively, and /£; and £, are the apparent activation 
energies corresponding to the initial A, and A, rates. 

If s represents the slope, (d log R)/[d(1/M,,)|, from 
eq (1 


Kk ! 


Kk, 9230.3. us 


Taking the common logarithm of this expression, and 
then the derivative with respect to the reciprocal of 
the absolute temperature, 


dlog K; 
d(1/T) aa 


dlog K, 
d/T) 


__@ log | <p (13) 


d(1/T) 


If each term is multiplied by —4.576107° (the 
molar gas constant appearing in the Arrhenius equa- 
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tion, 1.9872 X10™ keal mole~’ degree~', multiplied 
lhe —2.3026), eq. (13) becomes, for initial rates 
/ : ’ eal d log (—s) 

E,— E,= (4.576 X 107°) = , (14) | 


d(/T) 
| where E, and FE, are expressed in terms of kilocalories 
' per mole. ' 
From eq (11) one would expect k I 
stant within the experimental temperature range. 
Equation (14) could then be written as follows, for 
evaluation between temperatures, 7, and 7;: 


log : 
’ ’ ~— sd | ” 
iD E, 1.576 107”) ¥ (15) 
T, T. 
The apparent activation energies, /; and £), are 


based upon initial rates, and s should therefore be the 
initial slope. As remains constant at constant 


s 


to to be con- 


temperature, however, the accuracy of the value of | 


E.—E, calculated through the use of eq (15) should 
in the present case be greater than that which might 
be obtained by the usual method of determining each 
activation energy separately from the initial rates. 
The following experimental values of s are found 
for sample B: 13.79 < 104 at 300° C; —10.07X 10 
at 320° C;: and —9.04 10' at 340°C. The resulting 
values calculated for /4,— /, are 10.6 keal at 300° to 
20° C and 3.9 keal at 320° to 340°C. There seems 
to be no entirely satisfactory explanation for this 
difference, but the figure calculated for the 300° to 
20° C range may be preferred because the rates of 
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volatilization at 340° C are quite high from the ex- 
perimental standpoint. It cannot be ascertained 
from the present data whether a systematic experi- 
mental error is actually involved at 340° C or whether 
a change in mechanism is indicated. 

Cowley and Melville estimated the activation en- 
ergy for depropagation, E,, to be 18.5 kcal/mole from 
figures for the heat of polymerization (—13.0 kcal), 
and for the activation energy of propagation (5.5 
keal) [3]. If this value for #; and the value for 
E,— EF, in the range 300° to 320° C are substituted in 
eq (11), a value of 29 kcal/mole is obtained as an esti- 
mate of the activation energy for the transfer act in 
the degradation of sample B. 
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| Thermal Expansion of Binary Alkaline-Earth 
Borate Glasses 


Herman F. Shermer 


The thermal expansivity of a number of calcium, strontium, and barium borate glasses 


“ure reported 
ing alkaline-earth oxide concentration. 


minimum at about 20 mole percent of barium oxide. 
glasses were correlated with the density measurements in the liquid range. 


The expansivity of the calcium and strontium borates increases with increas- 
The expansivity of the barium borate glasses is a 


The expansivity measurements on the 
The expansivity 


of any composition is the lowest in the glassy (solid) region and is greater in the liquid region. 


The expansivity in a third region which was described as ‘“‘very viscous” 
was indicated to be greater than in either of the measureable regions. 


} but 


not measured 
The density 


was 


curves for each alkaline-earth borate series converges with increasing temperature. 


1. Introduction 


A study of the thermal expansion of the binary 
| akaline-earth borate glasses was undertaken as part 

of a comprehensive investigation of various physical 
properties of simple glass-forming systems. This 
report is the third of a series on the relationships 
between chemical composition and physical proper- 
ties of the binary alkaline-earth borates. The first 
report [1] ' was mainly concerned with the properties 
of the molten compositions, such as viscosity and 
surface tension, and the second [2] was concerned 
with the compressibility of the alkaline-earth borate 
glasses. 

The plan of the present work was the determina- 
tion of the relative effects of varving amounts of 
calcium, strontium, and barium oxides on the thermal 
expansion of borate glasses. It was desired to study 
these relationships over as wide a composition range 
as possible, but the region of glass formation was 
limited at low alkaline-earth oxide concentrations 
by the formation of two immiscible liquids, and at 
high concentrations by crystallization [4, 5]. 


2. Preparation of Glasses and Method of 
Test 


The glasses investigated were made by melting 
the required amounts of reagent-grade chemicals in 


platinum crucibles at approximately 1,300° C. 
The alkaline-earth oxides were introduced to the 
batch in the form of carbonates, and the boric 


oxide was added as boric acid. The melts, which 
weighed approximately 500 grams, were stirred with 
a propellor-type platinum stirrer until they appeared 
clear (at least 1 hour), and then they were poured 
into a cold iron mold. These binary glasses were 
later annealed by holding them at 590° C for 3 hours 
and then cooling them at 3 deg C per hour for 30 


Figur n brackets indicate the literature references at the end of this paper 
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hours, after which the power to the furnace was shut 
off, and the glasses were allowed to cool to room 
temperature. Polariscopic examination of the speci- 
mens revealed a slight amount of strain in the cal- 
cium borate glasses, but this strain was not detected 
in the expansion curves. The presence of strain is 
normally indicated by a decrease in the expansion 
with increase in the temperature in the range just 
below the point at which rapid expansion begins. 
The samples used in the determination of thermal 
expansion were cut from cylinders (approximately 
‘in. in diameter and 3 in. long) that had been used 
for compressibility specimens in previous work [2]. 

Bridgman and Simon [3] have shown that a com- 
paction of fused B,O;s occurs during compression 
measurements. As the glasses used in this study 
had been used previously for compressibility meas- 
urements and contained 63 to 84 mole percent of 
B.O;, one must consider the possibility that such 
compactions may have occurred. Weir and Shartsis 
[2] report a decrease in volume during compression 
testing on the order of 0.2 percent and also report it 
as being independent of B,O, concentration. An in- 
crease of this magnitude was attributed to experi- 
mental error. Even if this compaction were real, it 
is believed that its magnitude is too small to seriously 
affect the expansion data. Bridgman and Simon [3] 
have indicated that time and/or temperature may 
be employed to restore the specimens to their original 
densities. In the present investigation it is believed 
that the period between the compressibility measure- 
ments and expansion measurements (approximately 
6 months) and the heating of the specimens to 200° C 
preparatory to measurement minimized any possible 
effect of the previous compression measurements. 
An adjacent portion of the specimen was removed 
and analyzed by triple evaporation with HF and 
H,SO,." 


? The compositions were analyzed by Ernest B. Clark, formerly of the staff 
of the National Bureau of Standards. 








Following the original determination of thermal 
expansion, an adjacent portion of the four samples 
of calcium borate glasses that showed strain was 
reannealed to determine what differences, if any, 
would be observed in the thermal-expansion curves. 
These glasses were reannealed by heating them to 
655° C for 3 hours. They were then cooled at ap- 
proximately 3 deg C per hour until they reached 
455° C (68 hours), after which they were allowed to 
cool in the sealed furnace, with no further applica- 
tion of power. 

The measurements of thermal expansion were 
made by the interferometer method of Peters and 
Cragoe [6] as modified by Saunders {7}. The test 
specimen was formed into a T-shaped spacer about 
5 mm high, which was placed between two optically 
flat Vycor disks, the lower one of which served as 
an interference thermometer. This “sandwich” was 
placed in a metal cup with a fused-silica top to de- 
crease any temperature gradients. This cup was 
placed within a resistance-type electric furnace, 
which was heated at 2 deg C per minute from room 
temperature to above the deformation point. The 
motion of the two sets of fringes, one set formed 
within the lower plate, which indicated the temper- 
ature, and the other formed between the two plates, 
which indicated the expansion, was recorded on 
photographic film, and the results were read from 
this film. 


3. Results and Discussion 


3.1. B,O; 


In the study of the expansivity of the binary 
alkaline-earth borate glasses, it seemed desirable 
to know the coefficient of linear expansion and other 
characteristics of the expansion curve for boric oxide 
(B,O;) glass. Published data on the expansion of 
this glass are conflicting, as shown in the reviews of 
the subject by Morey [8] and Karkhanavala [9]. 
The differences seem to be caused primarily by 
different melting schedules and heat treatments and, 
as a result of this, differences in the amount of asso- 
ciated water [10]. Some of this water is retained 
from the boric acid, and other water may be adsorbed 
from the atmosphere after the glass is cooled. In 
order to reduce this water to an amount comparable 
with the other glasses measured in this investiga- 
tion, the boric oxide glass was melted on the same 
schedule as the binary glasses, and all glasses were 
stored in kerosine. Each sample was cut and dried 
by heating to 200° C and then allowed to cool in 
the furnace overnight, and the expansion was deter- 
mined the following day. The results of the deter- 
minations are shown graphically in figure 1 and are 
included for comparison in table 1. The small hump 
in the thermal expansion curve for BO; at about 
160° C is indicative of the small amount of water 
present [9]. 
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TABLE 1.—Thermal expansivity of binary alkaline-earth borate glasses 
| 
Chemical composition Density Critical Defor Average coefficients of expansion per degree C at—* 
; room come mation j j 
RO RO — te™P —_ Room to 100°C | 100° to 200° C | 200° to 300° C | 300° to 400° C | 400° to 500° C | 500° to 600° C 
Fused ByOy, 
Volk wt m ( ( 
0.0 0.0 1. 859 260 305 16. 1610 13. 77 X10- 
Ca0-B20, 
‘ =O { 184 633 658 5. 6410-6 6. OR & 10-¢ 6. 54x10 7. 1610" 7. 69X10-° 8. 68x 10-% 
’ - i Rea iled 646 662 5 49 6. 02 6. 39 7.13 7.71 8.09 
” O f 2. 531 647 665 5.73 6. 51 6. 74 7. 67 7.92 8.85 
1.4 . | Reannealed 649 665 5. 82 6 6. 81 7. 65 8. 10 8,31 
‘ “9 j 2. 5BS OAS 668 6. 01 6. 76 7. 52 7.78 8, 38 9.18 
4.9 _— Re iled Olt 6S 6. 62 6. 75 7 47 7. 85 9. X 9. 21 
. j 2 6 4 669 6. 40 7.04 7. 46 8, 22 &. 69 9. 55 
“ Re ile ‘ 669 6. 31 6. 84 7. 58 8. 16 8. 76 9. 35 
srO-B.O 
21.4 29.0 649 Us 614 79x10 6. 22 10- 6. 66X10- 7. 2810 &. 04 10-% 
21.8 2») 4 677 604 §21 77 6, 36 #. 7¢ 7.40 7. 87 
26. 0 44 2 202 19 635 RS 6. 46 7.03 7.76 8. 68 
. P { a5? 62 634s 6.01 6. 5S 6. 86 7. 56 8.12 9.04 10-6 
and LRe ired 621 638 97 6. 56 7.12 7.67 8.15 9. 08 
28 7 2. 931 624 41 6. 34 fi. 7¢ 7. 24 7. 98 8. 42 9. 42 
0) uF OOT7 2S 44 71 6. 96 ey 8, 28 9. 26 9.03 
42. ¢ 10 29 649 6, 96 7.47 8. 32 8. 49 8. 62 10. 16 
BaO-B,O 
17 Lf 2 G78 “wo 540 6. 5910 6. 69 10- 7. 25X< 10-¢ 7. 80x 10-¢ 9. 77X<X 10-6 
HK) 58 » I 1 S80 6.19 6. 59 7.01 7. 67 8.39 
z Mi). ¢ 2. 996 75 595 6.43 6. 64 02 7.71 8. 21 
27 5.1 oa S8u HOt 6, 59 6, 90 6.01 8.14 9. 36 | 
» { 17 ISU SS7 H06 6, 82 7. 46 7.89 9, 48 10. 39 
] w.6 $29 508 “10 7. 37 7. 64 8.13 8. 86 0. 58 
5.4 O4 H0o 7. 86 Ss. 48 & Ol 9. 6S 10.74 
» The estimated accura ! f average coefficient of expansion per degree C is +0.1X10 Ihe values given contain one additional figure for convenience 
intercomparisons and better determination of greater temperature ranges 
3.2. Calcium Borates calcium borates are included in table 1. Values for 
the reannealed samples are also given. In general, 


The results of the thermal-expansion measure- 
ments of the calcium borate glasses are given in table 
1 and figure 1. Although the four curves for the 
calcium borate glasses are similar, certain trends are 
apparent. The most obvious of these is the increase 
in expansivity with increasing concentration of 
calcium oxide. This is true in the region from room 
temperature to the ‘“‘eritical’”’ temperature and also 
at the deformation point. The critical temperature 
is defined here as the approximate temperature at 
which rapid expansion begins. This is in contrast 
with the expansion of fused B,O;, which may be 
thought of as having the lowest concentration of 
calcium oxide, and which has the highest expansion. 
It may also be noted that the temperature of the de- 
formation point increases with increasing calcium- 
oxide concentration. The deformation point is the 
highest point on an expansion curve obtained by the 
interferometer method and corresponds to the tem- 
perature at which viscous flow exactly counteracts 
thermal expansion. 

The values of the coefficients of expansion of the 
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the reannealed glass had a slightly lower expansion 
below the critical temperature. It may also be 
noted that the height of the curve at the deformation 
point was greater in the reannealed sample, indicat- 
ing a higher density at room temperature. The sam- 
ples used for reannealing were not large enough for 
accurate density determinations. 

The apparent increase in density of the calcium 
borate glasses when reannealed at 655° C indicated 
that the original annealing at 590° C was not con- 
tinued for sufficient time for equilibrium conditions 
to be established at that lower temperature. It also 
indicates that the equilibrium condition of the glass 
before and after the original annealing corresponded 
to a temperature above the equilibrium condition 
that was established during the second annealing. 
The results of the original measurement of thermal 
expansion are shown in figure 1. The results of the 
second determination are in table 1 but because of 
the small differences in expansivity found upon 
reannealing the curves are not shown. 
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Figure 2 shows the thermal-expansion curves for 
B.O; and the binary calcium borates plotted as 
change in density with temperature. ‘This method 
of presenting the data made possible the inclusion of 
density values in the liquid range that were obtained 
previously {1 It has not been possible to obtain 
density values of the liquids when they are very 
viscous, but it is hoped that the interpolations are 
reasonable. Three characteristics of this 
curves are the convergence of the curves with in- 
crease in temperature, the greater change in densits 
with temperature in the liquid range than in the 
solid range, and a middle range in which the glasses 
have higher coefficients of expansion than in either 
of the measurable ranges. 


set of 


3.3. Strontium Borates 


The thermal expansivitv values for 7 strontium 
borate glasses are given in table 1 and the curves 
of 4 of them are shown in figure 3. Most of the 
characteristics of the caletum borate curves are re- 
peated here. ‘The expansivity from room tempera- 
ture to the critical temperature and the expansivity 
and temperature at the deformation point all in- 
creased with increasing concentration of strontium 
oxide. ‘Two determinations of the glass containing 
26.9 mole percent of SrO are included in the table to 
illustrate the reproducibility of the results. 

Figure 4 shows the thermal-expansion curves of 
the strontium borates plotted as part of the density- 
versus-temperature curves for these compositions. 
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This set of curves shows a behavior similar to that 
of the caleium borates in its covergence with inereas- 
ing temperature, higher expansivity in the liquid 
range, and still higher expansivitv in the very 
viscous range, which was not measured. 


3.4. Barium Borates 


The results of the measurements of thermal ex- 
pansion of the barium borates are included in table 
1, and the curves of four of them are shown in figure 
5. There is one important variation in these results 
from those found in the calcium and strontium sys- 
tems. This change is in the relative expansion of 
the glass containing 20.1 mole percent of BaO below 
that of the glass containing 17.3 mole percent of 
BaQ. This decrease in expansivity with increase 
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q vids as a function of temperature 

in alkaline-earth oxide concentration is in the same 
relative direction as the decrease in expansivity 
between fused B,O; and the barium borate glass 
containing the least amount of BaOQ. The regions 
of glass formation of the ealcium and = strontium 
borates do not contain minima because of their small 
compositional range. The glass-forming region of 
the barium borates, however, is longer and contains 
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Figure 5 


the minimum. This minimum is better seen in 
figure 7. Although the expansivity of the barium 
borate glasses, below the critical temperature, passed 
through a minimum in respect to composition, the 
total expansion from room temperature to the defor- 
mation point did not. It may be noted here that 
the minima in the values of the coefficients of expan- 
sion with concentration do not occur at the same 
point in the various temperature ranges. In general, 
there seems to be a translation of the minimum to 
higher barium-oxide concentrations with increase in 
temperature. 

Figure 6 shows the relationship between densities, 
or expansivities, of the binary barium borates with 
temperature. The characteristics of this set of 
curves are similar to those discussed previously for 
the calcium and strontium borates. 
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3.5. Comparison of Thermal Expansivity of Alkaline- 
Earth Borates 


A comparison of the mean coefficient of expansion 
over the range from room temperature to 200° C of 
the binary alkaline-earth borate glasses as a function 
of concentration of alkaline-earth oxide in mole per- 
cent is shown in figure 7. It is apparent that a min- 
imum would have to occur in each of the curves if 
glasses could be made with lower concentrations of 
alkaline-earth oxides. There is not sufficient infor- 
mation available to determine the location of the 
minima in the caletum and strontium curves. 

Figure S&S shows a comparison of the deformation 
temperatures of the binary borate glasses as a fune- 
tion of alkaline-earth oxide concentration in mole 
percent. These curves show that at comparable 
compositions the calcium borate glass has the highest 
temperature at its deformation point, followed by 
the strontium and then the barium composition. 
This is in the same order as the melting points of the 
three alkaline-earth oxides. 
pa 


> 
4.'Calculation of Thermal Expansion From 
Composition 


Much work has been done in the past in an effort 
to find a series of factors from which thermal expan- 
sion could be predicted from the composition of vari- 
Sun and Silverman [11] compiled the 
results of several of these investigations. Among 
the oxides included in their compilation are three 


ous glasses. 


78 














700;-—_-_— ay 
anne 4 
ae 

o - 

*. 600 ai Bo 

Ww 7 : 

= y 

= 

dq 

a 4 

Ww 

a 

= 

Ww 

- 500 

2 

9 

- 

a 

= 

ec 

2) 

w 

Ww 

© 400 

« 8203 
300 
'@) 20 30 40 
RO, mole% 

FIGURE 8. Defor mation temperature as a function of COM posi. 
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BOs. 

\, calcium borates; 2, strontium borates; and CO, barium borat« 
that are applicable to the present investigation 


namely, BoO,, CaO, and BaO. The values assigned 
to these oxides by various investigators leave much 
to be desired. Figure 7 shows that the expansion- 
versus-composition curves of the binary alkaline- 
earth borate glasses are not linear, and, therefore, it 
would seem that in this case the coefficient of thermal 
expansion is not an additive property. 


5. Coordination of B.O, 


The thermal expansivities of alkaline-earth borat 
glasses are of special interest because they exhibit the 
so-called “boron anomaly’; 1. e., they show minima 
with respect to composition. Warren and Pincus 
[12] offer an explanation for this anomaly. It con 
cerns the ability of the boron atom to change from 
triangular to tetrahedral coordination when an 
oxide such as CaO is present and supplies the neces- 
sary oxygen, It is believed that fused BO, is tr- 
angularly coordinated, and that with the addition of 
CaO, increasing amounts of boron change to tetra- 
hedral coordination. 
demonstrated that bevond a point determined by 
the particular oxides present, any additional content 
of alkali or alkaline-earth oxides will tend to reverse 
the process, 1. e., change the BO, from fourfold to 
threefold coordination. A mathematical approach 
to the problem such as that of Warren and Pincus 
requires more X-ray data on the particular systems 
of glasses under study than are now available. 
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alternate approach is from the measurements of 
physical properties. The measurement of thermal 
expansion seems to be quite useful in this respect. 
The ‘ninima in the thermal-expansion-composition 


eurves shown in figure 7 are interpreted as indicating 
maxiina in the tetrahedral coordination of B.¢ ).. 
This seems to be in accord with the fact that the BO, 


tetraliedra tend to strengthen the structure and 
therefore have a lower expansion. 
6. Summary and Conclusions 
The thermal expansivity of glasses formed in the 


three binary alkaline-earth borate series was meas- 
ured. The thermal expansivity of the glasses in- 
creased with increasing concentration of calei ‘jum or 
strontium oxide. The expansivity of the binary 
barium borate glasses went through a minimum with 
changing concentration of barium oxide. A > mini- 
mum was implied in the other curves by the high 
relative expansion of fused BoO;. It was not found 
possiple to obtain additive fae tors for the calculation 
of thermal expansion that would be applicable to 
these compositions, 

The results were correlated with the density values 
obtained previously in the liquid range in order that 
the relationship between expansion and temperature 
could be studied over a greater temperature range. 
The expansivity of the solid glasses was the lowest, 
and that of the measurable liquids was much higher. 
The c ‘hange J Ihh de NSILV across the very Vv ise Ous region 
was such as to show a still greater expansivity in 
The density curves were found to con- 
In temperature, 


this range. 


verge with increase 


79 


7. References 


[1] Leo Shartsis and H. F. Shermer, Surface tension, density, 
viscosity, and electrical resistivity of molten binary 
alkaline-earth borates, J. Am. Ceram. Soc. 37, 544 
(1954). 

{2} Charles E. Weir and Leo Shartsis, Compressibility of 
binary alkaline-earth borate glasses, J. Am. Ceram 
Soc. (submitted for publication). 

[3] P. W. Bridgman and I. Simon, Effects of verv high 
pressures on glass, J. Appl. Phys. 24, 405 (1953). 

[4] P. Flint and L. 8. Wells, The svstem lime-borie oxide- 
silica, J. Research NBS 17, 745 (1936) RP941 

[5] Ernest M. Levin and George M. Ugrinic, The svstem 
barium oxide-boric oxide-silica, J. Research NBS 51, 
37 (1953) RP2430. 

16] C. G. Peters and C. H. Cragoe, Measurements on the 
thermal dilation of glass at high temperatures, BS 


Sci, Pap 16, 449 (1920) S393 

{7] James B. Saunders, An apparatus for photographing 
interference phenomena, J. Research NBS 385, 157 
(1945) RP1668 

[8S] G. W. Morey, Properties of glass, 2d ed., p. 273 Am 
Chem. Soc. Monograph No. 124 (Reinhold Publishing 
Corp., New York, N. Y., 1954). 

[9] M. D. Karkhanavala, Bibliography of thermal expansion 
of glasses, Glass Ind. 33, 404 (1952). 

10} James J. Donoghue and Donald Hubbard, Thermal ex- 
pansion studies of borie oxide glass and crystalline 
boric oxide, J. Research NBS 27, 371 (1941) RP1425 

[11] Kuan-Han-Sun and Alexander Silverman, Additive 
factors for calculating the coefficient of thermal ex- 
pansion of glass from its composition, Glass Ind. 22, 
114, 125 (1941 

[12] B. k. Warren and A. G. Pineus, Atomic consideration of 
immiscibility in glass svstems, J. Am. Ceram. Soc. 23, 
SOL (1940). 

[13] W. Stegmaier and <A. Dietzel, Die Bedeutung der 


von Glasschmelzen und Versuche zu deren 
Teil II, Glastech. Ber. 18, 353 (1940 


Basizitit 
Messung. 


WASHINGTON, November 1, 1955. 











oe 


ra 
ca 
ca 
0c 
sic 
th 
de 
sit 


th 
len 
rat 
tra 
gel 





' 


latt 
side 
the 
As | 
are 

(i 
on ¢ 

(I 
Int 

(¢€ 
to a 
It sl 
do I 
exhi 
beca 
dise 
plici 
tion 


1} 














11 of Research of the National Bureau of Standards 


Dielectric Relaxation for a Three-Dimensional Rotator in 
a Crystalline Field: Theory for a General Six-Site Model 
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2, February 1956 


Benjamin M. Axilrod 


A theory of dielectric relaxation is presented for a generalized six-site model where the 
transition probabilities for the turning of dipoles from one orientational site to another are 


arbitrary. The sites are arranged in three 
relaxation times generated by 
that the relaxation times are 
established 


1. Introduction 


Previous work [1]' with two simple models of a 
rather restricted nature has been sufficient to indi- 
cate that a set of discrete dielectric relaxation times 
ean arise in a lattice consisting of dipoles that may 
occupy orientational sites arranged in three dimen- 


sions. As in studies on the single-axis rotator 
(2, 3], the source of the set of relaxation times was 


the anisotropy in the crystalline field, which hin- 
dered the turning of dipoles from one orientational 
site to another. 

In the present paper the purpose is to extend 
the treatment of the three-dimensional rotator prob- 
lem «discussed in reference [1] by constructing the 
rate equations for a six-site model with arbitrary 
transition probabilities in order to examine the 
general nature of the set of relaxation times. 


2. General Six-Site Model 
2.1. Model 


A model in which a dipole of moment « on each 
lattice point ¢ ‘an orient in any six directions is con- 
sidered. These equilibrium orientations are ~_ 
the Cartesian coordinate 1 axes, sv, y, and 2 (fig. 1 (a)). 
As in previous work [1, 2, 3], the followi ing conditions 
are assumed: 

(a) There is no correlation in the motion of dipoles 
on different lattice points. 

(b) The barrier system is the same for all molecules 
in the lattice. 

The only elementary process is a reorientation 
to an adjacent site (single jump hypothesis, S-J/). 
It should be noted that the general results obtained 
do not depend upon the orthogonal site orientations 
exhibited in figure 1 (a). This model was selected 
because it serves as well as any other for the general 
discussion and it was desirable to derive certain ex- 
plicit expressions useful for polarizability calcula- 
tions for this special case. 


vures in brackets indicate the literature references at the end of this paper 


dimensions. 
this model is examined in a 
all real and positive 


| 
| 


| 
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The nature of the set of discrete 
general manner. It is found 


and the range of the relaxation times is 


The transition probability for the turning of 
dipole from a site 7 to an adjacent site 7 is given by 


— Wi 

k,j;=Be *7 (1) 

where B is a frequency factor, k is Boltzmann’s 
constant, and 7’ is the absolute temperature. The 


quantity W, is the local free energy of activation 
barrier, which must be surmounted by a dipole in 
turning from site 2 to site j (fig. 1 (b)). The transi- 
tion probability for the jump 7 to 7, k;,, is given by 
eq (1) with W,, in place of W,,. As transitions are 
permitted from a site 7 to four adjacent sites, there 
24 elementary process transition probabilities 


are 
k.; with 7, 7=1, 2, 6, and 7) in the general 
model. From physical considerations it is shown in 


the appendix (section 3.1) that of these 24 constants 
only 17 are independent, and that relations exist 
for the k,, corresponding to a closed circuit of tran- 
sitions and the inverse circuit. Thus, if we consider 
the series of site-to-site transitions 1-2, 2-3, 3-1, 
and the inverse circuit 1-3, 3->2, 2 


2, 21, we find ” 
Kyokogks, ky3kgoko1. 


Other relations such as the above must apply, 


(2a) 


one 


for each site triplet for the remaining seven octants: 
kisksaka =hiskasks (2b) 
ki gkaskss —hisksakar (2c) 
ksh sak) =hykosks: ; (2d) 


the four remaining equations, (2e) through (2h) are 
obtained from the equations above by substituting 
the subscript 6 for Any one of the set of eight 
relations in eq (2) may be derived from the other 
seven.’ These are the basic equations which re- 
strict the arbitrary choice of the &,; for the six-site 
model. Equations (2) are also needed for a mathe- 
matical proof noted later. 


2A single relation of this type was given in reference 3 for single-axis rotators 


see eq (7) and (25)). 
} For example, eq (2e), keokoskss=Kesksekes, may be obtained as the quotient [eq 
2a) X (2b) X (2e) X (2d) )/[eq (2f) x (2g) X (2h)]. 
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1, 2, and 6 Wie is the barrier for a transition from site 1 to site 2and WU © the 
barrier for the reverse transition t; and V2 are the free energies for sites 1 and 
2, respectively, relative to an arbitrary zero 


2.2. Rate Equations 


Assume a lattice with a total of NV dipoles where 
the number of dipoles oriented in the jth site at any 
time ¢ is given by N,. The summation of the N, 
is always "N. The number of dipoles turning per 
unit time from a given site 7 to an adjace nt site jis 


kN; and the number entering site + from j is kj,N 
the — leaving site 7 to turn to adjacent sites Is 
DkiN. Accordingly, the differential equations for 





; 


the change in population of each site are 


dN > 
“ hiethisthiuthisMi4 
ko Not+kaNe+kaNitkaN, 
= kieN, — (key + keea+Khes-+ lene) N 
, koNs+khoN,+keoN 
acid eee a ae ee a 
‘ kaN,+keN, 
, r (3) 
st eee a ee ee ae 
é‘ kN,+keN, 
= kN, +kasNo+kuN, 
(Kes +-Kenot+- Kes, -+- ens) Ng +e N 
sei kN. +kyeNs+keNi4 
Isp Ns— (kao t hos + oat hess) No | 





The molecular relaxation times, 73, will be included 


in the solutions of the above equations. 


2.3. Characteristic Determinant and Its Properties 


Solutions of linear homogeneous differential equa- 
tions of the type shown in eq (3) are of the form 


N= DiC pe", (4) 
8 


where the symbol /;(4) represents a function of some 
or all of the & The functions e”" are referred 
to as decay functions, since the f3(4) are always 
negative for B22 (see below). The index 8 refers 
to the mode of decay. Each decay funetion indi- 
cates the rate at which a particular mode of dlecay 
promotes the attainment of equilibrium after the 
abrupt removal of a disturbance. The f3(/) are 
obtained by forming the characteristic determinant 
of eq (3) and solving for its eigenvalues. This 
characteristic determinant is 


| j j he j 0 
h 1 hes 0 j h 
h h | ke 0) j 
AIDE, 
kev 0 ks, A, hee, 
j h 0) k,. A h 
0 h — kes, k, Bus | 
(5a) 
+ By the 23 terest is, / kos : 0 
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where we let D=d/dt and define A,, A», As, . . . as 
A; D- | t keyg his T hy; 
A D t ke, T kos hy; T hey, 
(5b) 
As D- ks, T kes T hes, T kee, 
Ae D t kee T Kees T kee, T kes. 


The characteristic determinantal equation, obtained 
by setting A equal to zero in eq (5a), has six roots, or 
eigenvalues, fs, for the operator D. The nature of 
the elements in A permits one to answer some 
questions regarding the eigenvalues and, hence, 
about the relaxation times, 73. 

The following questions may be raised: Are the 
eigenvalues real, what is their sign, what is their 
range, and are they distinet for the most general case, 
i.e., When the &,,; are distinct? The eigenvalues can 
be shown to be real by applying a theorem due to 
Goldberg.’ The line of reasoning is briefly thus: 
From the theory of matrices the characteristic 
determinant of a symmetric matrix has real roots. 
Hence, if a symmetric matrix, S, can be constructed 
that has the same characteristic determinantal 
equation as A, i.e., A=0, then the eigenvalues of A 
will be real. Sufficient conditions for forming such 
a matrix S with the same roots as A are the set of 
relations given by eq (2) and the inequalities 
k,>0. Since eq (2) and k,,>0 are imposed by 
physical considerations, it follows that the eigen- 
values are always real. 

The sign and range of the eigenvalues f3; are 
specified by applying a theorem derived by Ger- 
schgorin [4] and rediscovered by Brauer [5] regarding 
the areas in the complex plane in which the roots of 
a matrix lie; it can be shown (see appendix, section 
3.3) that for our particular case, eq (5a), the eigen- 
values lie between zero and twice the most negative 
diagonal element in the matrix (a,,) derived from 
Ain eq (Sa) by setting D=0 and changing the sign. 
Since, by definition, &,,>0, the roots are all negative. 

The proof that the eigenvalues are distinct if the 
k., are distinet has not been accomplished. 

From the theory of matrices we obtain two other 
properties of the characteristic determinant A. In 
the matrix (a,,;), defined above, it is evident (see eq 
(5b)) that each column sum vanishes. This in- 
dicates that one eigenvalue of A corresponds to the 
root D=0. Next, the trace or diagonal sum of a 
matrix (b,;) is equal to the sum of the eigenvalues, 
Dr,, of the characteristic determinant |),;—)é), 
where 6: is the Kronecker delta. From eq (5) it is 
seen that for the matrix (—a,,) the diagonal sum is 
>'k.,, where the prime means i#j. Since the eigen- 


values of the characteristic determinant are = un- 


———’ 





See appendix, section 3.2. 
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changed if each element is replaced by its negative 
the sum of the eigenvalues of A is —33’k The 


° J . 
latter property is useful as a check when the eigen- 
values are determined in a particular problem. 


2.4. Relaxation Times 


J 


With the general results obtained in section 2.3, 
we may write the solutions of the rate equations, 
noting that f; is zero, in the form 


(4’) 


where the zs are defined as the molecular relaxation 
times |1,2,3] 
ta= —1/fa(k). (6) 
Since we have shown that the fs are all real and 
negative, it follows that the 73 for the general six-site 
Cartesian axes model are all real and positive. As 
implied earlier, eq (4’) show how and at what rate 
a system tends to return to equilibrium after the 
abrupt removal of a disturbance such as a static 
electric field. Thus, for any arbitrary choice of the 
k,,, consistent with eq (2), the return of the system 
to equilibrium is described by an exponential decay 
and cannot include any oscillatory modes. 
Also, from the application of Gerschgorin’s Theo- 
rem to the range of the eigenvalues, we find that the 
shortest relaxation time Tmin is ° 


l 


mm 


> ; (7) 
2\a 


where @,, is the most negative diagonal element in 
the matrix (a,;) defined in section 2.3. If ky, is the 
largest of the &,,;, eq (5) and (7) lead to 


Tmin = gp. (3) 


pe 


The above results cover the inferences that may 
be readily drawn regarding the molecular relaxation 
times for the general six-site model. 


2.5. Dielectric Relaxation Spectrum 


If we assume the molecular and macroscopic 
relaxation times are equal, an approximation used 
previously [1,2,3,5], the molecular relaxation times 73 
describe the approach to equilibrium of the polariza- 
tion. The coefficients Cis (for B22) in general 
measure the “strength’’ of each mode contributing 
to the relaxation process. In the case where we 
consider the restoration of equilibrium following the 
removal of a static electric field, the polarization 


® In eq (7) the symbol |aam! represents the absolute value of the dam 








P(t) decays with time as 
PW =DP yo =D Pym, 9) 
r) 5 
The coefficients Ps, and hence the orientational 
polarizabilities as associated with —_ mode, mav 
be obtained from the quantities (822) for a 


particular model by standard ser Moe fl]. The ag 

are calculated assuming a polverystalline substance. 
We note that for sinusoidal fields eq (9) leads to 

the approximate expression [1,6] 

(10) 


€''(w) & KS agwrs/(1+w*r}), 


8 


where ¢’’(w) is the dielectric-loss factor as a function 
of angular frequency », A is a constant, and ag is 
the orientational polarizability associated with the 
8th mode. Although the expression is not exact, it is 
sufficient for an approximate analysis of the general 
characteristics of the dielectric relaxation spectrum. 
For the case of a single relaxation time, eq (10) 
reduces to the form of the Debye equation [7] for 
«’’(w). 

The coefficients Cjz (and hence the P3 and as) 
cannot be obtained exolicitly in the general case 
they depend (eq (3) and (4)) on the eigenvalues, and 
the latter cannot be expressed explic itly in terms of 
the k,,. As a consequence, we have not been able to 
prove in the general case, where the /,; are distinct, 
that all relaxation modes are active in the dielectric 
spectrum. However, the latter situation seems 
quite possible, since in earlier calculations inactive 
modes (compensatory reorientation of dipoles to 
maintain a net dipole moment of zero) were formed 
only for very simple models, such as those containing 
two pairs of equivalent sites. 

To obtain further informetion about the dielectric 
behavior of dipolar systems represented by the gen- 
eral six-site model, it is necessary to examine special 
models. This is planned for a subsequent paper. 

A few useful formulas to apply to special models 
are given in the appendix, section 3.4. 

The author is indebted to John D. Hoffman for 
numerous helpful discussions regarding the general 
subject and to Karl Goldberg and Martin Greenspan 
for assistance with mathematical aspects. 
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3. Appendix 
3.1. Relations Between k,,; 


Let the free energies relative to an arbitrary zero 
of energy for two adjac ent sites 7 and 7 be designated 
as V; and V,, and the free-e ‘nergy barriers for 2 to j 
and j to 7 transitions be defined as W,, and W,,, 
respectively. The height of the free-energy barrier 
between sites i and J; refe rred to the cero of energy, 
is then 

WO=Vit+W,=V,+U;, (lla) 


Consider a series of transitions 7 to /, j to /,/ tod for 
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a closed circuit of three adjacent sites and write the 
equations corresponding to eq (lla) for the barrier 
heights referred to the zero of energy 


jel WGYU)D=V;4+W,,=V.i4-W, (1 1b) 

lai WL) = Vt Wa =Vit+- Win. (11e) 
The sum of the eq (11) is 

Wt ,.4-W.=W,.4-0,4+ Wi. 12) 


If we recall the definition of &,;, eq (1), inspection 
shows that eq (12) corresponds to log (k,k)A),)= 
log (ki Ayky.), and hence 

kkk 


k gk yk ji. (13) 


The same type of relations, eq (12) and (13), can be 
obtained for any closed circuit of transitions between 
adjacent sites with three or more jumps to the cir- 
cuit.’ For the six-site model there are seven inde- 
pendent relations, such eq (13), namely, any 
seven of the eight eq (2). All others and those in- 
volving more sites are derivable from eq (13). Ae- 
cordingly, of the 24 transition probabilities only 17 
are independent. This is plausible on physical 
grounds. The local free-energy system of the general 
model is fixed when six minus one site energies and 
12 barriers Ws, 7) are specified; the “minus one” 
arises because the zero of energy is arbitrary. 

The following geometrical argument illustrates 
how all product relations of the type in eq (13), but 
involving more sites, can be constructed from suit- 
able products of the product relations given in eq (2) 
for three adjacent sites. 

The six orientational sites (fig. 1) are located at the 
corners of an octahedron. Associated with the 
product relations on the ky, eq (13), for a closed eir- 
cuit of transitions between the adjacent sites 7, j, /, 
is the cireuit of the perimeter of the face with vertices 
i, J, /, and the inverse circuit. Consider a closed 
circuit along the edges of the octahedron connecting 
the successively adjacent sites 7, 7, /, . for 
example, where there are 4, 5, or 6 sites. T his c ional 
circuit divides the octahedron, a simply connected 
surface, into two parts, each of which is constituted 
of triangular faces; for convenience, choose the sur- 
face with the fewer faces, labeled 7), . T,. If we 
traverse in succession in the same sense * as that on 
the whole perimeter (7,7, /,. . ., m, 7), the perimeters 
of each of the triangular faces f,. ., T, this is the 
same as traversing the bounding ‘path FS §, 


as 


7 An alternate method of arriving at product relations of the type in eq (13) or 
eq (2) is to invoke “detailed balancing”’ for each pair of adjacent sites. This 
means that at equilibrium the number of dipoles turning in unit time from site 
i to site j is equal to the number turning from site j to site i, or kiyjCa=kyCa; Ca is 
the equilibrium number of dipoles in the ith site. For a particular case, sites 1, 2, 


, we may then write kywCy =k Cn, f Cus kon, and ky kioCny. The product 

of the last three e¢ quations yields eq (2 a), Kisktovh koh 
§ That is, if the perimeter of the partion of the surface bounded by the sites 
.. m, i) is traversed in a clockwise sense by the above sequence, then each 


Sede constituent triangles is traversed in the same sense 
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m, i. Hence if we multiply the product equations 
such as (13) for the triangles 7;,. . ., 7,, traversing 
their eiges in the correct sense, we obtain the corre- 
sponding product relation for the cireuit (7, j, J, id 
m, i) and the inverse circuit (7, m, "<4 i). ‘In 
this wits the key produc t relation for any sequence of 
connected sites can be formed by multiplying the 
product relations for the enclosed triangular circuits. 


3.2. Application of Goldberg’s Theorem * to the 


Matrix for the Six-Site Model 


Theorem: If a real matrix A= (a,,) of order n has 
the properties a;,2>0 for 1,j7=1,..., n with 14) 
and 

4 i,4isi, - - - Ui, =U,i,4i,i, ~~ - Vii, (14) 
for all (7 , im) belonging to(1,...,n),m>3,then 


all the characteristic roots of A are real. 

We wish to apply the theorem to the matrix — Apo 
derived from A in eq (5a) by letting D be zero and 
changing the signs of the elements. If we let a,; be 
the element of — Ao in the 7th row and the jth column, 
the usual order of index notation, we note that for 
ix} 


Qy=k;, (15) 


the order of indices is reversed on the k’s. Henee 
eq (14) can be written as follows, after 


1. €., 
with eq (15), 


interchanging the left and right sides and then re- 
versing the factor order on the right side: 
hi, Vivi ‘ hi,i; ki,iskixi, os kisi, 
htsinn ‘ hisi,ki.i,. (14’ 


Equation (14’) now evidently equates the product of 
ky’s for the closed circuit of m (3<m<6) 
hn ho, lo— 1g lL» hy to the corresponding 
product for the inverse circuit i: hens <a do- } 
There are two types of circuits: (a) circuits in which 
at least one successive pair ip, te are nonadjacent sites, 
and (b) circuits in which all successive site pairs are 
adjacent. Since by the S-J hypothesis ie 
ki; =0 for all nonadjacent site pairs, eq (14’), and 
hence eq (14), is valid for case (a). Next we have 
indicated in section 3.1 that closed-circuit product 
relations of the type eq (14’) are valid for a sequence 
of three or more adjacent sites. Hence our matrix 
— Ay for the pose six-site model satisfies the above 
theorem. 


sites 


act 


3.3. Range of Eigenvalues 


Brauer [5] has derived 
regarding the characteristic 
matrix, 


the following theorem 
roots (eigenvalues) of a 


*K. Goldberg of the Applied Mathematics Division of the Bureau. His proof 
$ published in RP2652 in this Journal. The theorem as stated above is in a 
slightly different form from that as stated by Goldberg; his condition (2) has 
deen replaced by the condition aj;>0 w hich still pe rmits the construction of an 


equivalent symmetric matrix, namely, ig —(a4j054) 2. 
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jth circle (7=1, 


Theorem: Let A=(a;,) be an arbitrary matrix and 


n 
p> a,\/=P;, dSlar,|=Q,). 
= i=1 
jut ixj 


Each characteristic root \ lies in at least one of the 
circles 


<P; 
and in at least one of the circles 
|z—aj,| <Q. 


A is a square matrix of order n over the complex 
field; z is the complex variable x+-.y. P,is the sum 
of the absolute values of the off-ciagonal elements in 
the ith row; Q; is a similar sum for the jth column. 
In the matrix —Ap derived from the characteristic 
determinant, eq (5a), by letting D be zero and chang- 
ing the signs, we know that the off-diagonal elements 
are all real and positive, and so an off-diagonal 
column sum is Q;. Also from eq (5a) we see that 
the jth diagonal element aj; is —Q;. Hence, for the 
, 6) the radius is Q,, and the 


center is on the real axis at —Q,;. Let Q,>Q, 
(7=1, , 6). Q, defines the largest circle. Pre- 


viously all roots were shown to be real. Hence all 
roots lie on or within the mth circle on the axis of 
reals, that is between zero and 2a,,,, (Anm being the 
most negative of the diagonal elements in — Ao). 


3.4. Useful Formulas 


Relation between the equilibrium number of dipoles 
in each site with and without a static electric field. We 
refer to the method of solving the rate equations in 
reference 1. After a “mode analysis” in which for 
each given mode 8 the coefficients Cis are determined 
relative to one another, the boundary condition at 
t=0 is introduced. This condition corresponds to 
the sudden removal of a static electric field. This 
results in the equations 


i ‘i> — 2 "8 


where ("is the equilibrium number of dipoles in the 


ith site in the presence of the field. A general 
expression for these coefficients in terms of the 


equilibrium coefficients for zero field is helpful in 
solving any particular model. 

The equations below for the C4 are derived in a 
manner similar to that indicated in section 2.4 of 
reference 1. 


- 

We assume a static electric field F defined by the 
magnitude F and the polar angles & and @, with & 
being the azimuth angle measured from the z-axis. 
For convenience, let 


y=uF/kT and y, 


Yu 


y sin &— cos @ 
y sin — sin @ 


Yz=Y Cos &. 
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A Matrix with Real Characteristic Roots 
Karl Goldberg 


It is proved that a certain matrix, which is the coefficient matrix of a differential equa- 


| tion found in the theory of dielectric relaxation, has only real characteristic roots. 


This 


is done by finding a real symmetric matrix with the same principal minors and thus the 


same characteristic roots. 


In papers [1, 2, 3]* by J. D. Hoffman and B. M. 
Axilrod a certain differential equation has a real, 
constant coefficient matrix A=(a,,) with the prop- 
erties 
for all 7 


ad >( 


iy SO for all i] (2) 


all 7 


a ay 


for all ,.,;., (4) 


In this note we prove the conjecture of J. D. 
Hoffman and B. ML. Axilrod that a matrix with 
these properties has only real characteristic roots and 
that these characteristic roots lie between zero and 
twice the maximum diagonal element. 

The second statement is a direct consequence of 
the first and properties (1), (2), and (3) by the 
theorem of S. Gershgorin [4] and A. Brauer [5]. 
This theorem states that all the characteristic roots 
of a matrix cl=— (a lie in the area bounded by the 
a < +> a a 





circles 


Therefore, we may concentrate on proving that all 
the characteristic roots are real. We shall do this by 
exhibiting a svmmetric matrix with the same char- 
acteristic roots as one with properties (2) and (4). 

Let B= (b be a matrix of the same order as the 
matrix sl (a which has properties (2) and (4), 
such that 


4 (d,.da u/e 


with the square root assumed to be positive. Then 


1b h 


—— 


' 


Now consider any principal minor 


<1 (81,82, ij, i,j Sm of A. 


motes 
Krom the definition of a determinant we have 


* 
A(8;,82, Sm) S CAs 2,4s,t, ° 


t m 


e 9 


Sim 
where the sum is taken over all permutations 
S So 
-m ...4 


m 


and the sign is positive if the per- 
mutation is even, and negative if it is odd. 


Of 8,82, . . . . Sms 


Each permutation is the product of eveles. Thus 
| ; 
SiSo . 1 e Sm / l i i, 
- ? 
tit, .. bol dy 
so that 
a 4 aty =m ai, ofl aig Ci, iy 
h. bd , Cas 
1 ig” izitg ip ty 
b..1,6 » Gidtudlie m'm* 
It follows that 
AA(8),82 . . . 8m B(s,,8, . . . 8m), 


that is, that the corresponding principal minors of A 
and B are equal. 

This implies that 4 and B have the same charac- 
teristic equation and thus have the same character- 
istic roots. 


{1} J. D. Hoffman and B. M. Axilrod, Dielectric relaxation 
for spherical molecules in a crystalline field: Theory for 
two simple models, J. Research NBS 64, 375 (1955 
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[5] A. Brauer, Limits for the characteristic roots of a matrix, 
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Reflection and Transmission of Gamma Radiation by 
Barriers: Semianalytic Monte Carlo Calculation ' 


Martin J. Berger and John Doggett 


The transport equation for photons which have been Compton-scattered repeatedly with 
a specified sequence of energy losses and deflections was solved analytically, by an ‘‘orders of 
scattering”’ approach, to obtain conditional reflection and trans‘nission probabilities. These 
probabilities were then averaged over all possible intermediate angular and energy configura- 
tions in a Monte Carlo calculation carried out on the NBS automatic computer (SEAC), 
The efficiency of the Monte Carlo calculation was increased by taking advantage of the close 
correlation between the transmission of radiation through thick barriers and the correspond- 
ing penetration problem in an infinite homogeneous medium for which an exact solution was 
available. The reflected and transmitted energy flux (integrated over all spectral energies 
and directions) has been calculated for 0.66-, 1- and 4-Mev radiation incident on water barriers, 
and 1-, 4-, and 10-Mev radiation incident on iron, tin, and lead barriers. Angular and 
energy spectra were obtained for water at 0.66 Mev. The ratio of the scattered energy flux, 
transmitted through a barrier to the scattered flux at equal depth in an infinite medium, was 
found to be constant for distances from the source greater than four mean free paths of the 





= 


source radiation. 


1. Introduction 


The theory of the diffusion and penetration of 
gamma radiation in an infinite homogeneous medium 
is now well developed, complete ‘solutions bei ing 
provided by the moment-method [1] * and the ¢ losely 
related semiasymptotic Laplace-transform method 
(2). A satisfactory theory of boundary effects is still 
lacking. The “orders of scattering’ approach [3] 
provides formal solutions to boundary problems, 
which can in practice be evaluated only when few 
orders of scattering need be considered. The Monte 
Carlo method [4] is a very effective means for treating 
boundary effects but is inefficient when applied to the 
radiation that penetrates very far from the source. 
This inefficiency can, in principle, be removed by 
so-called “biased”? sampling, but the success of such 
techniques has not yet been demonstrated con- 
clusively. 

This investigation is part of an exploratory Monte 
Carlo program for the survey of boundary effects in 
gamma-ray diffusion. The specific aim is to caleu- 
late the reflection and transmission of radiation by 
plane parallel barriers. The approach is eclectic, 
and a combination of all three of the methods de- 
scribed above is being used. It can be characterized 
as an “orders of scattering’? method in which the 
multiple integrals, which ordinarily make the evalua- 
tion of the solution troublesome, are evaluated 
numerically by random sampling. Conversely, one 
could eall it a “reinforced” Monte Carlo calculation 
whose efficiency for dealing with deep penetration 
has been increased in two ways: by confining the 
random sampling to the angular and energy variables, 
while treating the space variable analytically, and 
through the use of correlated sampling, the cor- 
related problem being the diffusion of y-radiation in 
an infinite homogeneous medium, 


This work was supported by the Office of Naval Research and the Reactor 
Divis sion of the U. S. Atomic Energy Commission. 
? Figures in brackets indicate the literature references at the end of this paper 


The reflection and transmission problems differ 
substantially from the point of view of random 
sampling, because transmission through thick bar- 
riers, being unlikely, tends to require an excessive 
amount of sampling and thus puts a high premium 
on analytical help. Nevertheless, the two problems 
were treated together because they fitted in a single 
computational scheme. 

Two types of questions may be asked in radiation 
diffusion problems. One concerns the estimation of 
the radiation intensity at any one position, cor- 
responding to the observable response of some 
integral detector like an ionization chamber. The 
other type of question is more detailed and concerns 
the spectral composition and directional distribution 
of the radiation at any one point. In the work 
reported here most of the effort has been devoted to 
intensity (energy flux) determinations. Exploratory 
studies of spectral and angular distributions were 
made only for one source energy (0.66 Mev) and one 
barrier material (water). 

For the deep-penetration problem the attention 
was focused on the comparison of the radiation 
intensity behind a finite barrier of thickness ¢ and 
inside an infinite medium at a depth z (equal to @) 
from the source plane (see fig. 1,a,b). No quantita- 
tive indication was previously available on this 
comparison. Notice that practical situations are 
usually intermediate between those schematized in 
figure 1,a,b, because there is some matter (possibly 
of different composition) behind the barrier, though 
often at some distance from it (fig. 1,e). 

On the problem of relative intensity a rather 
definitive answer has been obtained. In the first 
place, it proved unnecessary to have extended the 
calculation to very great depths because the in- 
tensity ratio no longer varies appreciably beyond 
about four mean free paths of the incident radiation. 
It follows that any future investigation may confine 
itself to moderate depths and may thus utilize less 
elaborate methods. Data on the values of the 
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FIGURE 1 Schematization of the barrier reflection-transmission 


proble in 


intensity ratio for various depths are given in table 3 
for four materials (water, iron, tin, and lead) and 
four y-ray energies (0.66, 1, 4, and 10 Mev Those 
data may be applied as correction factors to the 
extensive results of calculations for infinite media 
[5]. The ratios differ substantially from unity only 
for low energies and low-Z materials, and also vary 
smoothly, so that adaptation to other circumstances 
ought to be easy. 

In the treatment of the reflection problem, atten- 
tion was directed in the first place to the dependence 


of the reflected intensity on the y-ray energy, on 
the material of the barrier, and also on the angle of 
incidence. The results for a semi-infinite barrier 
are given in table 2,a. An additional item of 
interest was the dependence of reflection on the 
barrier thickness, that is, on how rapidly the re- 
flection approaches the value for a semi-infinite 
ee, Some data on this aspect are given in 
table 2, 
2. Outline of the Calculation 
2.1. Analytical Frame 

We consider a homogeneous plane-parallel layer 
of material located between the planes :—0 and 
-=/. Radiation is incident with specified energy 
and direction on the face z=0. What are the 


characteristics of the radiation reflected through the 


plane 0 and transmitted through the plane {? 
The state of a photon traversing the barrier is 
deseribed by three variables: the energy Kk the 


| 
| 


90 





the 


angle @ between the direction of motion an 
2-axis, and the depth 2 in the barrier. 

Photons can be scattered and absorbed inside 
the barrier. Let yuy(/) equal the probability of 
absorption per unit path length; ¢(/2) the probability 
of scattering per unit path length, and ¥(/.6; /’ 6’ 
the probability distribution (normalized to unity) 
that a scattering will change the state of the photop 
from (42,6) to (E’ 0’). Consider a photon that, as 


the result of m scatterings, has successively assumed 
the energies and directions 
Bok, « «+ ley nergy-angle history / 
energyv-angie history /,. 
Mins « oe on ; 


The index n refers to the state of the photon imme. 
diately after its nth scattering; »=0 pertains to the 


state in which the photon enters the barrier. The 
probability distribution of h,, is 
" I ,. , 

Q,= TI WE, On; Em+19m+1) (1) 


Let P,,(2,t4,) be the probability that a photon with 
history A, will cross the plane 2 (O<2<¢) after the 
nth and prior to an (n+1)st seattering. An an 
alvtical expression for this crossing-probability js 
derived in section 3. If a beam of A’ photons per 


unit area and unit time is incident on the barrier. 
the flux per unit time at depth 2 is 
F(E,0,z,t)=K Sof, | ae 6, : dE, 


| 10,2°.Qu8 k—E,)s(0—0 (2) 
a 
where 6 is Dirac’s delta function. If the weight 
factors f, are unity, the flux represents the number 


of photons crossing a unit area of the 2-plane per 
unit time. Most radiation detectors do not measure 
the number flux through a fixed plane; their response 
depends in various ways on the energy and direction 
of the radiation and is represented by the /,’s. A 





few examples are: 


yy Fl : 
“lux r ocurrent Tr unit 
a | Number Flux (or current) per uni 
: area of surface perpen- 
b E, Knergs ; 
dicular to the z-axis 
Flux per unit area of sur- 


face perpendicular to the 


direction of motion of the i 
c sec 6, Number radiation This is the 
d sec 6, E, hnerg,y flux seen bv an Isotropic | 
detector It is propor- ' 
tional to the photor 


density 


The numerical results of this paper pertain to flux } 
tvpe (al). By this choice they become comparable to 
the results of an extensive of caleulations for 
infinite media by the moment-method [5] 
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2.2. Monte Carlo Estimates 


In the evaluation of expression (2) for the flux, the 
sum over collisons must be cut off for n->N, and the 
integration over the intermediate energies and angles is 
earried out by the Monte Carlo method. <A suitable 
value of the cutoff N depends on the absorbing prop- 
erties of the medium, the depth of penetration, and 
the spectral-energy range over which one wants to 
know the flux. Calculations for an infinite medium 
show that for all practical purposes the energy flux 
yanishes below 30 kev in water, below 50 kev in 
iron, 79 kev in tin, and 100 kev in lead, for source 
energies between 0.5 and 10 Mev. The photon 
energy histories were therefore terminated at these 
limits. Trial calculations indicated that under these 
circumstances, for a barrier thickness up to 16 mean 
free paths (uof=16), it was sufficient to let N=25 
for water, 12 for iron, 9 for tin, and 6 for lead. 

Energy-angle histories of length N were sampled 
according to the distribution function Qy. The 
detailed equations for doing this are well known. 
They are presented elsewhere, together with the 
adaptations required or desirable for computations 
on an automatic computer [6]. 

Let h,,; be an initial portion of length n<.N of the 
jth sample history hy,, and let P,,,(2,t) be the corre- 
sponding crossing-probability. If a sample of J 
histories is available, the flux of radiation in a 
specified spectral interval 2, and angular interval 
R, can be estimated as follows: 


. » a. N 
| JE {dP EO, 2.1) =+ >>» 
Re J Rs ‘ ln=0 
Singh nj (2b) Ar (En3) A223), (3) 
where 
AE f—1,if E,. isin the interval R, 
oe a if fe, is outside the interval R, 
(0 f=1.if @,,isin the interval I, 


{=90, if 6,, is outside the interval 2s. 


If one wants to know the value of the flux for fixed 
values of FE and @ rather than in intervals, the 
sampling procedure must be modified so as to be 
conditional not only on the initial photon state but 
also on the final state #,—F and 6,=—86. For this 
purpose one replaces P,, by P*(AEW(E,_;,9,-1; FE, 9), 
where 


We have been content to obtain the flux in energy 
and angular intervals because the main body of 
calculations was done in any case for the total 
intensity integrated over all spectral energies and 
directions. 

In conformity with the usage that has become 
common in shielding work, the results for the 
intensity are presented in the form of buildup 
factors (>t), which are equal to the ratio of the 


total to the unscattered energy flux: 


E "© 
| dE | d0F(E.6.2.t) 1) 


KfoPo 


The Monte Carlo estimates of the buildup factors 
are the mean values 


B(z,t)=(1 JD) SSBy(2t), (5) 
j=1 
where 
- 
B ;(z,t)=(Eo| see 00| Po)~' 55 E,;| sec 6,5! P,;(2,0). 
n=0 
The standard deviation of the buildup factor 


estimate is J~'/*e(z,t) where 


1 J -_ 1/2 
a(2,t) s 3 18,(2,0—B(=or }- (6) 
J j=l 


The buildup factors (and their standard devia- 
tions) were estimated from the same set of energy- 
angle histories for three specific situations: 

(i) B(O,t) energy reflection buildup factor 
(ii) B(t,t) energy transmission buildup factor 
(iii) B(t,©)=energy buildup factor for the flux at 
depth ¢ within a semi-infinite 
medium. 


3 


2.3. Correlation 


The difference between Bit) and B(t, @) is pre- 


_ dominantly due to photons with histories such that 


in an infinite medium they would cross the plane 
t, have their directions reversed as the result of 
one or a few collisions in the vicinity of t, and would 
then make another crossing. Thus S(t,) is 
loosely speaking—a function of two factors: a prob- 
ability p,(t) of penetration to depth ¢, and a proba- 
bility p. of being turned around in the vicinity of the 
t-plane. On the other hand, /4(¢,t) only contains the 
factor p,(t). The greatest statistical fluctuations in 
a Monte Carlo calculation are found in the estimate 
of p,(t), because only a small part of the sampled 
histories provides information concerning deep pene- 
tration. If the same set of energy-angle histories is 
used for the estimate of B(t, ©) and (t,t), it is to be 
expected that the estimate of the ratio of these two 
quantities will have a much smaller statistical dis- 
persion than the estimates of the individual buildup 
factors because the ratio no longer contains the 
troublesome factor p(t). Unscattered radiation 
should be excluded from the comparison because it 





The energy reflection buildup factor has the’ ollowing physica! significance. 
Consider an isotropic detector placed in the path of a broad beam of radiation, 
ind denote the energy flow through the detector by Fo. Then place a barrier of 
thickness tf immediately behind the detector. This increases the energy flow 
through the detector to F. Then #(0,1)=F/ Fo. The quantity B(0./)—1is related 
to but not identical with the energy albedo, which as commonly defined pertains 
to an energy current. 








i 

s not affected by the boundaries. We therefore 
subtract 1 from both buildup factors and introduce 
the ratio 


Bit} 


. (7) 
Bit, o)—1 


k(t) 


estimate of k has an estimated 


26. where 


i+ 


A Monte Carlo 


standard deviation J 


ie _ olb,t) 
Bit.t) 


a(t,@) 


+ |- 
L Bit, @) ] 


a(t,?) al(t,a@) 1/2 
2p — — ’ (S) 
Bit,t)—1 Bit, ©)—1 

and 

J a a 

1/J) >> (By (t,) —Bit,d] |B (t, ©)— Bit, @)] 
p(t) — (9) 
a(t, t)a(t, o) 

is the correlation coefficient of A(tt) and Bt, @). 


The two buildup factors have a strong positive cor- 
relation so that p~1 and the value of o,/k is small. 
For moderate and deep penetrations, /)(f,©) is 
practically identical with the buildup factor for an 
infinite medium which can be obtained by a reliable 
moment-calculation. These known results, together 
with an estimate of the ratio /(t), provide an ac- 
curate indirect method for determining the trans- 
mission buildup factor which requires much smaller 
sample sizes than a direct estimate. Thus the 
Monte Carlo and moment-methods complement 
each other: one provides knowledge of the boundary 
effect, the other knowledge of deep penetration of 
radiation. In section 5 the advantage of this cor- 
relation technique is demonstrated. 


3. Solution of the Transport Equation for the 
Crossing-Probabilities 


The crossing-probabilities obey the following 


integral equations: 


P,.,(2) e(E,) dP, (2 )\em Wau) +k 
' cos 4, J0 
(10a) 
if cos 6, Ly 0. 
clE,) "¢ s 4 P 
P. (2 — dz’P.(2’) [wal Ens +c} ji< COS One 
cos @, P 
(10b) 
if cos 6,.,<0. 
Pf z)= eal Bo +e (100) 


These equations express the fact that the probability 
of passage through the 2-plane after n+-1 collisions 
is equal to the probability of passage through the 
plane 2’ after n collisions, times the probability of 
a collision between 2’ and 2’+<dz2’, times the probabil- 
ity of passage from 2’ to z without further scattering 





| or absorption, the triple product being averaged 
over all values of 2’ compatible with the boundary 
conditions. , 

We shall find it convenient to introduce in place 
of the probabilities P, the related quantities 


G,(z)=|sec 0,|P,,(2) (11) 
which already incorporate the trigonometrical factor 
required for the type of flux which we want to deter. 
mine. The following abbreviated notation will be 
adopted: 


uy, —alE,)+e(E,) ) 





c¢,=E,) 
N,—= Cos 6, "a (12) 
¥.=G if cos 6, >0 
u,—=tif cos @,<0 J 
| Equations (10a, b, and ¢) can then be rewritten in 


the form 


+ oe & FE Sate tty 
G,, «\( o) " ™ d- Gz Me ta (13a) 


"n+l 
(13b) 
No 


To solve these equations, assume that G@, can be 
expressed in the form 


(14) 


y Cy nN 4 ad Nin 
G, (2) - 2 
Nn+1 m=0 Nn+1K Nm B 
\ “= , wm Me " rk Uh } (15) 
ié f ft \ 


so that the functional form of G,,, 
that of G,. The expansion coefficients are 


( l . , 
= At forO<m<n (16a) 
n+l +1 My 
Nn+1 7 
{n+} _S n+1 ~ (161 
“eo 7% — As ( ” yD) 
m ti 


We have thus obtained inductively the general 
solution for G,(z). Although the solution is rather 
formal and not readily surveyed, it is convenient for 
numerical evaluation with an automatic computer 
because of its recursive nature. 


92 








is the same as | 


1 
Jari 
But 
bee: 
inte 
pall 
by 3 
of | 
bot 
equ 
At 


If t 
ced 
For 
con 
limi 
fori 
intr 


ax 
by 
evel 
eXal 
10 


1 
aul 
two 
of ¢ 
sam 
The 
sec’ 
to 
Any 
obt: 
(r,s 
outy 
1 
tion 
n 
Kit] 
ean 
stal 
as t 
in t 
\ 
the 
of 
dist 
T 











12) 


‘nD in 


13a) | 


13b) 


(14) 


tion 





(15) | 


ie 7 


| 6a) 


6b) 


eral 
her 

for 
iter 


It is possible that mptins1— Mn4idm=O, so that singu- 
Jarities occur in the recursion relations (16a, b). 
But this difficulty is apparent rather than real, 
because upon substitution of the coefficients Az*! 


into the expression for G,,, the singularities occur 
pairwise and cancel each other. This can be seen 
by a reexamination of (15). Note that the existence 
of the singularity implies that 11%: >0, because 
hoth uw. and wp) are positive. It follows that w,, is 
equal (0 Ungi- Therefore, the factor multiplying 
An in eq (15) is 


C,¢ ~ 
a= é "= | (17) 
Mn+l ’ 
n+l . 
AS Mmbn41— 7 M 0) 
—_” m 
a”"—c,(2—u,)e ™ : (18) 


If the singularities are removed by a limiting pro- 
cedure, the functional form of G,,, is thus changed. 
For the numerical computation on an automatic 
computer we have found it simplest not to use the 
limiting procedure, but to stick to the unmodified 
formal solution, while eliminating singular cases by 
introducing a cutoff for the quantity a@= (sys)! p11) 

Un’ Mm) In the denominator of eq (16a). Every 
a<10-° was replaced by 10 The error incurred 
by this approximation increases with |2—u,, But 


even for deep penetration it is quite small. For 
example, when a(2—4,,)= 20, [exp(— 0.001 20)— 1]/ 
10 20.201, so that the error is only 1 percent. 
4. Machine Computation 
The calculation was programed for the NBS 


automatic computer (SEAC). It was carried out in 
two stages. First, energy-angle histories (in groups 
of approximately 100) were generated by random 
sampling and stored on magnetic wire or tape. 
These histories were then used as input data for the 
second stage, in which the quantities G@, (proportional 
to the crossing-probabilities P,) were calculated. 
Any desired type of flux or buildup factor can be 
obtained by a suitable linear combination of the 
the selection being made by an appropriate 
output code. 


Gi. 3. 


The buildup factor code provides for the computa- 
tion of B04), Btt), and B(t,e), for poof =2"d 
n=O,1, , 5), 4 bei ing a disposable parameter. 


Either all 18 or any subgroup of these buildup factors 
ean be computed in an integrated operation. The 
standard peviiions a(0,t), o(t,, and a(t, ), as well 
as the correlation coefficient p (t) are also contained 
in the output. 

Most of the computational effort was devoted to 
the buildup faetors, but limited use was also made 
of output for energy spectra and angular 
distributions. 

The amount of SEAC computing time depends 
on the number NV of collisions through which one 


codes 
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follows the photon histories. As has been mentioned, 
it was found sufficient to let N=25 for water, 12 for 
iron, 9 for tin, and 6 for lead. The amount of 
machine time (seconds per history) required for the 
generation and storage of photon energy-angle his- 
tories was 
t,=0.4 N. (19) 
For the joint computation of 18 buildup factors, 
the machine time (seconds per history and buildup 
factor) 
4,=0.18 N+0.0205 N? (20) 
for 6<N<25. 
buildup factor, 
tory) 


For the computation of an individual 
the machine time (seconds per his- 

&=1.3 N. (21) 
It can be seen that the joint computation of several 
buildup factors was more economical. 

A major difficulty in setting up the machine pro- 
gram was the sealing problem. SEAC operates only 
on numbers less than 4 in absolute value. But the 
coefficients A® in the expansion for G, can occasion- 
ally become extremely large because of sampling 
accidents. This difficulty was aggravated by the 
fact that singularities were not removed analytically 
from the recursion relations (l6a, b) but by a 
numerical cutoff procedure. This resulted in the 
occurrence of pairs of coefficients with approximately 
the same large absolute value, but opposite sign, 
whose weighted difference must be calculated with 
great accuracy. This required the use of floating- 
decimal-point routines, which, although they were 
specially adapted to the problem, proved to be 
costly in computation time. 

The occurrence of large terms with alternating 
signs in the series for G, provided a convenient check 
on the code and the computer operation. Clearly 


G,, being a positive multiple of a probability, must 
be a positive quantity. It was found by experience 
that because of the alternating character of the 


series for G,, almost any sort of code or machine 
trouble quickly led to the occurrence of negative 
G,’s. The machine was instructed to interrupt the 
computation whenever this happened. # As many 
G,’s are computed per second, the absence of such 
interruptions for prolonged periods was a reasonable 
guarantee that the operation of thefcomputer was 
error free. 


5. Results 


5.1. Statistical Efficiency 


Before presenting the main body of our results, 
we shall discuss the statistical accuracy of the Monte 
Carlo estimates in detail for a sample calculation 
based on 300 photon histories, pertaining to 1-Mev 


radiation incident normally on an iron barrier. In 
table 1 are listed, in addition to the estimated 
buildup factors and buildup factor ratios, other 


quantities of statistical interest including fractional 
standard deviations and correlation coefficients. 














TABLE 1. Statistical analysis of buildup factor estimates for 1-Mev radiation incident on an iron barrie: TAR 
Jo 
ri 
: a(Q), t - a(l.t =, ; a(t,= ; j eK mi 
os Bio. ¢ Riot ; Btt Bit ; Bit, « Ru ; p ‘ ] 
0.5 1. OAS 0. OSS6 1. 40 0. 61 l. 53 0. 43 & 2X 1@-2 0. 755 0. 29 0.4 
1. 0 1. O57 . 0696 1. 71 73 1. 89 . 58 2.2 10°? . 798 . 20 3 
2.0 1. O61 0747 2. 43 1. 00 2. 68 89 6.5X 10-3 851 15 15.2 
1. 0 1. O61 . O7AT 1. OF 15S 1. 45 1.48 1 9X 10-3 . 890 IS 1. 6 V 
8. 0 1. O61 O7A7 7. 80 2. 90 8. 60 2. 68 3. 4X 10-3 895 31 67.2 
16. 0 1. 061 0747 17.8 187 20.0 58s. 85X10" 884 1 15.0} 
I; 
' 
Now let us see what the large amount of analysis | TaBiLe 2. Energy reflection | 
introduced into the Monte Carlo calculation has 
accomplished. To estimate the transmission buildup 1. Energy Reflection Buildup Factors for Semt-infinite Barriers L 
factor for a barrier with thickness pof=16 with a | a - 
standard deviation amounting to 5 percent of the Water | iron rin Lead F 
buildup factor would require on the order of 10° . 
photon histories in an ordinary random sampling | *"™ 
calculation. But according to table 1 the semi- - - . = ‘ = ‘ “i 
analytical method requires only 8.5 10° histories. 
Taking into account the estimated increase in Mer. 
computation time per history required in the semi- of AE 
analytical calculation, one finds that the effective 1.0 1.081 1153 1.061 1.142 -1.022) L086 1.009 1.042 
reduction factor is ~2,500. This is respectable but ms | tan) tae wen — 
not quite good enough for a computer with the 
capabilities of SEAC. But one can do better still b. Dependence of Reflection on Barrier Thickness * 
through the use of correlated sampling, 1. e., by ) . 
estimating the buildup factor ratio A(é). It can be mol 
seen from table 1 that, in accordance with expecta- wanes ees ° | Ir 
tions, the correlation coefficient p is very close to : : ' 
unity, and the fractional standard deviation ok "~ 
correspondingly small. A 5 pereent standard devia- Water 0. 66 0 0.65 0.88 0.99 : 
tion now only requires a sample of 67 histories. — on = — - ~~ 
ia hl . . . ° “ r 
faking into account the increased computation ae 1.0 o 4 os 1.0 I 
time per history, the reduction of effort is S80-fold. Lead 1.0 0 "7 1.00 1.00 
The over-all improvement factor is therefore 
y 000 SO = 200,000. » Columns 4, 5, and 6 give the values of k(0) =[BO,t)—1] [Bio, © 1] 7 
The last column of table 1 lists the factor q¢ pe | tetnyepren ttt! dp nr gh eae, Co 
loft O/BAOP/ ok, by which the correlated sam- tame 
pling technique reduced the required sample size. | the source radiation is already equivalent to a semi- | wou 
Statistical analyses at other energies and for other | infinite medium. sho 
scattering media yield over-all improvements in Table 3 contains an extensive list of buildup | bv : 
efficiency of about the same magnitude. factor ratios k(t), for water at 0.66, 1, and 4 Mev, for | mor 
iron, tin, and lead at 1, 4, and 10 Mev. Again, it? bes 
can be seen that the boundary effect, indicated by than 
5.2. Buildup Factors deviations of k(t) from unity, is appreciable only at} med 
low energies and for low-Z materials. It is signifi- | sma 
; ‘ . eet cant that the ratios stay constant, within the limits | toq 
Reflection buildup factors for semi-infinite media of statistical error, as the barrier thickness is in- | T 
are listed in table 2,a, for diverse conditions. The creased from of =4 to wot =16. The physical reason | the 
wee" sizes were such that the quantities B(0,t)—1, | for this constane v is the establishme nt of radiative Mol 
i. ¢., the increase of the flux due to reflection, had a | near-equilibria in regions far from the source, which | by t 
sti atistic al accuracy of 5 percent. It can be seen | is well known from “the diffusion theory for infinite , thie 





that energy reflection is appreciable only for low | media. Sufficiently deep in the material, both the ing 
energies and low-Z materials. Table 2,b, shows the energy spectrum and angular distribution of the | set 

depe ‘ndence of reflection on the barrier thickness for | radiation flux become very slowly varying functions | infir 
some typical cases. The outstanding feature of these | of the penetration depth. Under these cireumstances | witl 
results is the rapidity with which the reflection of | it is not surprising that the boundary effeet, which | tran 
energy assumes its maximum value. A barrier with | depends on the characteristics of these distributions, | of 6 
a thickness of one, or at most two mean free paths of | likewise remains practically unchanged. , fors 
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TABLI Buildup factor ratio k(t) 
for « paring transmission of radiation through a finite bar- 
rier th penetration to an equal depth in a_ semi-infinite 
mea mnt 
pot 
Ma il Energy 
1.0 1.0 0 1.0 s.0 16.0 
Mer 
| 0. 66 0. 601 0. 663 0.713 0. 783 0. 785 0. 7§4 
Wat l 661 720 74 S21 SJ . S30 
| 4 .s4u ASS 912 G20 26 433 
| l 7) TUS S5l SUO) . 895 su 
Irot 4 SuO) 910 423 436 32 ay 
| 10 “41 59 . 972 u74 97S 977 
] . SSY Y1l y24 935 USS 46 
Tin { 441 426 YAS W67 74 Q7s 
119 951 60) Y62 973 971 969 
| ] a3y 951 Y“6Yy 975 _u79 . O82 
Lead 4 4) “77 Ys2 yuo) uur UO4 
| 10 ye G40 yyus uy U4 ~ 005 
Estin ited 
wecuracy 
0 >) 1.5 l , 0 2. 
* Comparison of a barrier with an infinite medium, 
TABLE 4. Energy-transmission buildup factors 
peat 
Energy 
Material Ver 
iT) 10 2.0 10 sO 16.0 
| 0. 66 1. 40 1. U6 3.10 wy 13.3 34.4 
W l 1. 40 1. 80) 2.72 Ol 10.5 25. 7 
l 4 122 1.42) LS) 260 $21 © 7.20 
| 1 1. 40 1.72 2.43 10 7.80 17.8 
I ‘ 1. 2 1. St 1.72 2 it 4.17 7.4 
| 10 107 11 1.35) 1.7 2. 80 s. 
| l 1.29 1. St 2.10 3.1 41 «10.2 
] ‘ 1. 16 1. 31 1. 3 2.4 4.12 9. 41 
| Ww 106 1.12 1. 26 1 ; }. &. 22 
| ! 1.20 1. 35 1. 63 2.09 2.8 41.24 
Lea i 111 1. 23 1.44 1. US 3. 28 7. 46 
Lx 1.03 | 1.08] 1.17] 1.40 2.17 | 647 


The results of table 3 pertain to a comparison of 
finite barriers with a semi-infinite medium. For dis- 
tances of >4, a comparison with an infinite medium 
would vield practically identical results. Column 3 
shows the buildup factor ratios for pof=1 obtained 
by a Monte Carlo calculation for transmission and a 
moment-caleulation for an jnfinite medium. It can 
be seen that these ratios are up to 10 percent lower 
than the corresponding ratios for a semi-infinite 
medium. For got =2 there is a similar but extremely 
small reduction, but our data are insufficient for us 
to quote a value of the reduction. 

Table 4 lists the transmission buildup factors for 
the same range of conditions as in table 3. For 
ut =0.5, 1, and 2, the results were obtained entirely 
by the Monte Carlo calculation. For greater barrier 
thicknesses the values quoted are the result of apply- 
ing the correction factors of table 3 to an extensive 
set of calculations for an infinite medium [5]. The 
infinite-medium results are stated to be accurate to 
within 5 pereent. The accuracy ot the derived 
transmission buildup factors is therefore of the order 
of 6 to 7 percent. The direct Monte Carlo results 
forshallow penetration have about the same accuracy. 
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5.3. Application of the Monte Carlo Method to an 
Infinite-Medium Problem 


For heavy elements the boundary effect is so small 
that the Monte Carlo method does not add much 
new knowledge. Yet it can be utilized as a con- 
venient and laborsaving tool for the extension of 
available results. As an example of such an appli- 
cation we have considered the penetration of radia- 
tion from a plane monodirectional oblique source in 
an infinite medium. This problem can be handled 
by the moment-method, but requires much more 
work than the corresponding problem for plane 
monodirectional sources with zero obliquity, or 
plane isotropic sources. We considered a 4-Mevy 
plane oblique source in lead. <A set of 100 energy- 
angle histories were selected that vielded buildup 
factors in close agreement with those of a moment- 
calculation for a monodirectional zero obliquity 
source for penetrations up to wot=16. Then we per- 
formed a rotation of the angular part of these his- 
tories. Provision for doing this is built into the 
SEAC code for the computation of the G,’s. After 
rotation, these histories can then be used as input 
data for oblique source problems, and the results can 
be expected to be reliable, since the energy part of 
the histories, known to be properly representative, 





= 
NORMAL 
INCIDENCE 








0000! 
ie) 


FIGURE 2. j- Mev plane 
monodirectional oblique Source through a le ad harrier. 


Transmission of radiation from a 








remained unchanged by the rotation. In this man- 
ner we have made quick calculations for obliquity 
angles 0)=45°, 60°, 75°, and 90°. As a final check, 
the fluxes from sources with different obliquities 
were summed so as to obtain the flux from a plane 
isotropic source. Again the agreement with a direct 
moment-calculation was close. Thus a fix on the 
accuracy of the Monte Carlo results was obtained at 
two points, and in spite of the small sample size the 
penetration of radiation to very great depths could 
be caleulated with an accuracy which we estimate 
to be 5 to 7 percent. The results of the caleula- 
tion are shown in figure 2 in the form of plots of 


the percent-energy-transmission (=buildup factor 
exp(—wpol)) versus pof, for different obliquity 


angles. 


5.4. Energy Spectra and Angular Distributions 
In figure 3 we show for normally incident 0.66- 
Mev radiation the energy spectra of the scattered 
flux transmitted through water barriers with thick- 
and yot=8. The corresponding energy 
spectra at the same depths in an infinite medium 


NESSES Lol 


are also shown. (The histograms pertain to the 
Monte Carlo results, the curves to a moment- 
calculation.) The shaded areas between corre- 


sponding histograms for transmission and an infinite 
medium are a measure of the boundary effect. It 
can be seen that the presence of boundaries removes 
a large part of the low-energy end of the spectrum, 
but leaves spectral components near the source 
energy unchanged. It is noteworthy that the 
spectra are quite similar at both depths. This lends 
support to our argument that the constancy of the 
buildup factor ratios is connected with the occurrence 
of radiative equilibria. 

Figure 4 pertains to the same conditions as figure 
3, and shows the angular distribution of the scattered 
transmitted energy flux, and of the corresponding 
flux in an infinite medium. The shaded area 
again indicate the boundary effect. The shaded areas 
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Ficure 3. 
throug! 
mally incident beams of O#6-Mev photons. 


Enerqy spect a of scalle red radiation transmitted 


water barriers of thickness pot , and pot —8, for nor- 


The unshaded parts of the histograms pertain to penetration through a finite 
barrier, While the shaded areas indicate the additional radiation that would pene 
trate to a depth ¢ in an infinitely extended medium rhe curves accompanying 
the histograms represent the results of calculations by the moment-method for 
an infinite medium, 





| for negative values of cos@ represent photons which 
| crossed the boundary 2, 4, 6, times, while the 
areas for positive values of cosd represent photons 
which crossed the boundary 3, 5, 7, - - - times in ap 
infinite medium. The jalatively small number 
of photons in the latter group indicates that repeated 
crossings of any plane deep in an infinite medium 
are quite unlikely, so that it should be possible to 
calculate the boundary effect in reasonable approx. 
imation by infinite-medium theory, by setting the 


transmitted flux equal to the infinite-medium | 
flux integrated over positive values of cosé. The 
angular distributions are similar to the chergy | 


spectra in that they do not change their shape ag 
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Figure 4. Angular distributions of scattered radiation trans 


mitted through 
for normally incident beams of O.66-Mev photons. 


water barriers of thickness pol=4 and yp g 


The unshaded parts of the histograms pertain to penetration through a finite 
barrier, While the shaded areas indicate the additional radiation that would 
penetrate to a depth ¢ in an infinitely extended medium. The angle @ measures 
the direction of the transmitted radiation with respect to the normal to the 
barrier face. Shaded areas for positive values of cos # represent photons that 
would cross the boundry plane at least three times in an infinite medium. 
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Figure 5. 
tion reflecte d from a semi-infinite 
incident beam of O.66-Mev photons, 


water barrier, for a normally 


rhe angle @ indicates the direction bet ween the direction of the reflected photons 


ind the normal to the barrier face. 


Enerqy spectrum and angular distribution of radia 








f 


the 


pol = 
EF 


tion: 
inci 
Bot! 
loca 
exist 
asso 
mat 
app! 


A 
and 
expe 
tion 
Mor 
resu 
of t 
boul 
feret 
thro 
mate 

Tl 
twee 
mea: 
valu 
fact 
resul 
so tl 
wish 

iN 
pene 
thro 
mad 
alun 
dens 
per 
as fe 
ries, 
med 
eula 
resu 
effec 
the 
smal 
the 
One 
to a 
eXPe 
ther 
ing ¢ 
in tl 
mue 
pass 
side) 
dept 
as ft 
fron 
cale 








vhich § 


© the 
otons 
In an 
mber 


pated 


dium § 


le to 
Prox. 
r the 
dium 

The 
lergy 


i 


1 finite 
would 
sures 





.0 


adia- | 


’ 
naily 


otons | 


ee 


as | 








LI 


the barrier thickness is increased from jof=4 to 
t=8. 

Figure 5 shows the angular and energy distribu- 
tions of reflected radiation due to 0.66-Mev photons 
incident normally on a semi-infinite water barrier. 
Both distributions have two distinct peaks. Their 
locations are such that they appear to imply the 
existence of two preferred modes of reflection, one 
associated with a first scattering through approxi- 
mately 180°, the other with a first scattering through 
approximately 90°. 


bey 


6. Discussion 
6.1. Comparison with Experiments 


A recent survey of experiments on the penetration 
and diffusion of y-radiation [5] indicates that few 
experiments have been done for plane source radia- 
tion that are comparable with our calculations. 
Moreover, the uncertainty of the experimental 
results is usually at least of the order of magnitude 
of the boundary effect. It would seem that the 
boundary effect could best be checked by a dif- 
ferential experiment in which one measures the flux 
through a detector, both with and without backing 
material behind the detector. 

The penetration of radiation with energies be- 
tween 0.66 and 2.76 Mev through iron slabs has been 
measured by Beach et al. {7].. Their experimental 
values agree to within a few percent with the buildup 
factors listed in table 4. But the experimental 
results include large corrections for air-scattering, 
so that the evidence is not as conclusive as one might 
wish. 

Kirn, Kennedy, and Wyckoff [8] measured the 
penetration of Co” radiation (1.17 and 1.33 Mev) 
through concrete barriers. An attempt has been 
made to match their results by a calculation for an 
aluminum barrier with a concrete-equivalent electron 
density. (The photoelectric absorption cross section 
per electron is very nearly the same for concrete 


as for aluminum, except at extremely low’ ener- 
cies.) The penetration in an infinite aluminum 
medium was found by interpolation from cal- 
culations of Goldstein and Wilkins [5]. These 
results were then corrected for the boundary 
effect with the use of table 3 of this paper. Finally, 


the predicted buildup) factors were increased by 
small amounts (not more than 2%) to account for 
the energy response of the experimental detector. 
One more major adjustment must be made prior 
to a comparison of theory and experiment. The 
experimental cobalt source was 1 em thick, so that 
there was an appreciable chance for Compton scatter- 
ing and resultant energy degradation of the radiation 
in the source. The degraded radiation is absorbed 
much more strongly than the primary radiation while 
passing through the concrete barrier. When con- 
sidering the attenuation at moderate and great 
depths one should renormalize the source strength so 
as to include only undegraded photons emerging 
from the seuree. An approximate single-scattering 
calculation leads to the estimate that the effective 
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source strength is .82+.05 times the nominal source 
strength. Table 5 shows the experimental and 
theoretical energy transmission (for normal incidence 
of the Co® radiation on the concrete barrier). The 
agreement with the adjusted experimental values is 
good. Although the source strength correction 
constitutes a weakness of the comzarison (at least 
insofar as verification of the boundary effect is con- 
cerned), it is significant that by adjusting a single 
normalization constant one can get agreement be- 
tween theory and experiment over a wide range of 
barrier thicknesses. 

TABLE 5. Energy transmission by a concrete barrier 


Comparison of the experimental transmission found by Kirn et al. [8], and the 
predicted transmission. (Transmission—=buildup factor Xexp(—pol).) 


Energy transmission, in percent, 
with concrete barrier thickness of 


5.0in. 10.0in. 15.0in 17.5 in 

Experiment, uncorrected 39.0 11.0 2. 80 1. 40 
Experiment, corrected for the deg- 

radation of radiation in the source 4.8 13.2 3. 36 1. 6S 

Theory 44.3 13.7 3. 41 1. 64 


The experimental literature on the reflection of 
gamma radiation from barriers is even scarcer than 
that on transmission. We have learned of some 
unpublished measurements by Kirn, Kennedy, and 
Wyckoff * of the reflection of Co” radiation from 
effectively semi-infinite concrete and lead barriers. 
(The lead barrier actually was a concrete block 
fronted by \ in. of lead, but this combination was 
in effect equivalent to a solid lead barrier.) By 
interpolation based on the results in table 2 we 
obtained reflection buildup factors which, according 
to the comparison in table 6, are in reasonably good 
agreement with the experimental values. 


Energy-reflection buildup factors for Co™ radiation 
incident on semi-infinite barriers 


TABLE 6 


Comparison of the theory with the experiment of Kirn et al 


Angle of a ’ a 
Material incidence, @ Experiment rheory 


dea 
. { 0 1. 063 1. O05 
Concrete | ww 1. OS6 1.110 
| 0 Lol 1. OOUS 
Lead | 60 1. 041 1. 040 


6.2 Comparison With Other Calculations 


In table 7 a comparison is made for buildup 
factors for shallow and moderate penetration depths, 
pertaining to 0.66 Mev radiation incident on water 
barriers. The comparison includes the results of 
this paper, infinite-medium calculations by the 
moment-method [10] and the results of two other 
Monte Carlo calculations, one of them done by a 


direct stochastic analog method [9], the other 
‘HH. O. Wyckoff, private communication 








ZABLE 7. 
methods for 0.66-Mev radiation in wate: 


Comparison of buildup factors obtained by various 


medium 


Reflection 
from semi- 


Penetration in 


Pransmission semi-infinite 
} v) ° 
Method Furs medium, — 
af, — medium, 
Bi, = 
B 
mM 
2 ‘ ! 2 ‘ 0 wo 
Monte Carlo this 
paper 1 $i 5. Wy 102 4. 1d 1. 148 1. 189 
NMionte Carlo (NBS 
[6] 1.89 3.21 (462.48 407 698 115 118 
Monte Carlo(NRL 
Y) 1.92 3.30 255 2.43 . 75 6. 55 1.14 
Spencer-Fano mo 
ment-method [10 389 7.37 


based on an analytical integration over a sampled 
collision density [6]. The agreement of the various 
methods is on the whole satisfactory. 


6.3. Comments 


We have demonstrated the constancy of the 
boundary effect for thick barriers, and the resultant 
lack of necessity to do deep-penetration Monte 
Carlo calculations, only for plane geometry. But it 


appears plausible that a similar situation will also 
prevail in other geometries since the constancy of 
the boundary effect uppears to be rooted in the 
occurrence of equilibrium spectra far from a radiation 
source. 

The advantage of incorporating knowledge ob- 
tained from other sources into a Monte Carlo caleu- 
lation through the technique of correlated sampling 
has been clearly demonstrated. In fact, it was only 
this technique which made it possible for us to survey 
the boundary effect for a variety of conditions with- 
out excessive computation. 

The analytical treatment of the space variable 
likewise proved successful, but we feel that it may 





| have been too much of a good thing. From the stand. 
point of computing efficiency, there is an optimum 
amount of analysis that should be put into a Monte 
Carlo calculation, which depends on the nature of 
the problem and the capabilities of the computing 
machine. It is questionable whether we operated 
near this optimum. If one interested in the 


Is 


radiation intensity, it is not unreasonable to take, } 


as we have done, relatively small samples and to 
analyze them exhaustively. But if energy spectra 
and angular distributions are desired, it is probably 
preferable to use larger and more representative 
samples of energy-angle histories and to analyze 
them less thoroughly. 
would be to confine the analytic treatment of the 
space variable to a limited number of collisions, and 
then to revert to a stochastic calculation. 


The authors thank Anne Futterman for aid and 
advice in programing the problem for SEAC. 
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An Examination of the Helium Vapor-Pressure Scale of 
Temperature Using a Magnetic Thermometer’ 


E. Ambler and R. P. Hudson 


The variation of the mutual inductance of two coils surrounding a paramagnetic crystal 
has been measured as a function of the saturation vapor pressure of helium in the range 1.3 
to 4.2° Kk. The fact that this quantity should vary inversely as the absolute temperature 
has been made use of to investigate the consistency of two recently proposed vapor-pressure 
temperature scales. The results suggest errors above 2° K in the empirical equation pro- 
posed by Clement, Logan, and Gaffney (in contrast to the experiences of Erickson and 
are in closer accord with the thermodynamic 


Roberts with a magnetic thermometer) and 


calculation of Van Dijk and Durieux. 


1. Introduction 


All practical thermometry in the “‘liquid-helium 
region’ (that is, between 1° and 5.2° KX) is dependent, 
directly or indirectly, upon a knowledge of the rela- 
tion between the saturation vapor pressure of helium 
and the absolute temperature. Many investigators 
measure directly the pressure over the liquid-helium 
ervostat and derive T from p-T tables; others use 
the same procedure to calibrate a resistance ther- 
mometer or magnetic thermometer, etc. In this type 
of vapor-pressure measurement the accepted practice 
is to apply a depth correction (“hydrostatic head” 
to the measured value of p; this correction is, will 
ever, of somewhat doubtful validity and can be 
avoided by measuring the pressure in a vapor-pres- 
sure bulb that is in thermal equilibrium with the 
material under investigation (effectively zero  im- 
mersion 

The p-T tables in general use at the present time 
comprise the “1948 seale” [1]? and are based on 
the work of Schmidt and Keesom [2], Bleaney and 
Simon [3], and Kamerlingh Onnes and Weber {4}. 
Possible errors in the 1948 seale were admitted at 
the time of its preparation, notably in the 1.83° to 
22° K region from consideration of the helium 
isotherm measurements of Kistemaket (5), and be- 
tween the normal boiling point and the critical point 
due to the sparseness of the experimental data upon 
which the seale was based in the latter region. The 
investigations of Erickson and Roberts [6] with a 
magnetic thermometer for the region 1° to 4.2° K 
and those of Berman and Swenson |7| above 4.2° kK 
with a gas thermometer provided a strong basis for 
arevision of the 1948 scale, and an empirical formula 
was developed by Clement and coworkers {8} from 
which a p-T table could be conveniently calculated 
to any desired precision, and which fitted the new 
data to within 0.002 deg throughout the entire range. 
Asummary of the situation obtaining in October 1954, 
prior to the development of the Clement formula and 
the publication of the supporting evidence obtained 
by Corak et al. [9] from calorimetric work and by 

A brief account of this work w presented at the Fourth International Con 


ference on Low Temperature Physic Paris, August 30 to September 8, 1955 
Figures in brackets indicate the literature references at the end of this paper 
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Keller [10! from He’ and He® isotherm data, has 
been given by Hudson [11]. The latter report made 
reference to some preliminary magnetic thermometer 
investigations by Hudson and de Klerk (unpublished ) 
that were in qualitative agreement with the findings 
of Erickson and Roberts (6| but suggested some what 
larger errors in the 1948 seale in the region of 3.5° K. 
These measurements have now been extende d. polos 
improved apparatus, and the results provide the 
subject of the present report. While this work was 
in progress, Keesom and Pearlman |12]| reported that 
some anomalies in their calorimetric data could be 
removed upon reevaluating their data in terms of the 
Clement equation. 

In summary, the evidence available in mid-1955 
strongly supported the validity of the Clement equa- 
tion. A revision of the 1948 scale, long overdue, 
therefore seemed feasible and the time opportune 
with the approach of the Fourth International Con- 
ference on Low Temperature Physics in Paris in 
September 1955. During this conference, however, 
the results of a new thermodynamic calculation of 
the p-7 relation, which differs by several millidegrees 
from the Clement equation below 1.5° K and above 
2.2° K, were announced by Van Dijk and Durieux.’ 

The first part of section 4 of this paper deals with 
an analysis of our data in terms of the Clement equa- 
tion. In section 4.2 the same data are reanalyzed 
in terms of the Van Dijk-Durieux table. 


2. Apparatus 


Po reduce uncertainties in the measurement of p, 
the saturation vapor pressure, occasioned by meas- 
urement of the bath-pressure and application of the 
hvadrost: atic-head correction, a vapor-pressure bulb 
was emploved. A series of measurements was per- 
formed with the apparatus shown in figure 1, a. 

ry © ° 

The magnetic thermometer comprised a paramag- 
netic salt specimen, A (a l-in. sphere ground from a 
large crystal of chromic methvlammonium alum)! 

H. van Dijk and M. Durieux, paper presented at the Fourth International 
Conference on Low Temperature Physics, Paris Aug. 30 to Sept. &. 1955. We 
ire indebted to Dr. van Dijk for kindly providing us with advance information 
on this ealeulation 

* The choice of salt was based upon the requirements Reproducible behavior 


mall erystal field) splittin ind dipole-dipole interaction, absence of direct 
exchange interaction, and ease of growth of larve erystals 








“e A, — 10-* deg; at 2° K, 1.2107 deg; g 


°K,5& 107° deg. This precision was not re «alized “a 

| in cae tic e because of small fluctuations at the lowes; ‘s 
| pressures and to the necessity for applying smal “ 
corrections for the nonlinearity of the bridge, by s 


the probable error due to such effects was less thay th 











al 1 millidegree thoughout the range of measuremen é 

17] | (1.3° to 4.2° K). | ob 

G. | IL] Subsequently the apparatus was rede signed, fof th 

Mi iil of reasons given below, and the second version is showy ap 

li} U4 in figure 1, b In this modification there was pp it 
yy | | vacuum jacket shielding the tube, F, the vapor} 4c 
sini pressure bulb was shortened conside rably, and access ms 

] Ale to the latter was made possible by winding the coi hy 

/ |) of the mutual inductance on demountable formers alc 
| B 1|~ This assembly proved thoroughly satisfactory from} g; 

yi’ a = the point of view of rigidity. ° 

al 

IZ, | i * 
Zin 3. Experimental Procedure = 
Ud} os {P With liquid helium in the cryostat, liquid helium} P° 
| i} } Seolod was introduced into the vapor-pressure bulb by a 
it D | | condensation of gas under a small overpressure | ji¢ 
ui} ett | Prior to the silvering and final assembly, a test filling (1- 
ih was carried out to determine the exact — ft no. 

|. {I helium required to submerse the crystal, A, and to co 
E | observe the change of liquid level within the bulh 
(a) C (b) during continuous cooling of the bath. This level _ 

will tend to rise because of condensation of the helium f 
4 FriGure 1, Vapor-pressure bulb and magnetic thermometer. gas in the manometer sv stem and fall bee use of pe 
., First apparatus, used in June 1955 measurements; b, modification as used in | Increased density of the liquid as the lambda point “we 
August 1955 measurements, For explanation of symbols, see text is approached. In the present apparatus, the former sa 
effect tended to outweigh the latter, and ther - 
resulted about a 3-mm submersion at the lambda ete 
located within a mutual inductance, P-S, which was | point. As the precise depth varied with bath levd[ ga 
wound, in a rigid assembly, on concentric glass tubes. | (a large part of the dead-space gas is contained in th in 
The latter was connected to an a-c mutual-inductance | section of F that is immersed), the possible correction pone 
bridge [13]. For convenience of removal and inspec- | of 0.03 mm of Hg was not applied to the measured} jy 
tion of the erystal, the bulb, B, was extended to | value of p. [The corresponding error in 7 is 3 10" apy 


project beyond the coil windings, and closed at the | deg at the lambda point and becomes progressively} mo 
bottom by a ground-joint plug, C. A pedestal, D, | smaller as 7 increases; below the lambda point ther pre 


supported the crystal, and the volume of the bulb | is, of course, no depth correction. ] r 
was reduced by means of a hollow glass ‘‘filler’’, E. In an experimental determination of helium vapor} yj 
The tube, F, connected the bulb to the manometer | pressure the most probable sources of error, i.e) gp, 
system and was protected by a vacuum jacket, G. | extraneous heat influxes to the bulb, will lead to an} gj, 


Sealed-in glass disks, each perforated by one small | overestimate of p. Apart from the more familiar} fp), 
hole (about 2 mm in diameter), provided radiation | radiation and conduction effects, there is one ~~ the 








shielding, and F and B were silvered internally. to low-temperature apparatus; viz., a hea influx the 

The purpose of the jacket, G, was to avoid a | due to oscillations in the gas column in a tube such} to 1 
“cold spot” on F: Due to the hydrostatic head effect, | as F. As a result of these heat leaks, the surface ‘ 
the liquid helium in the surrounding bath is coldest | temperature of the bulb liquid will be raised, and} 9}, 
at the surface and, if the vapor in F comes to thermal | temperature inhomogeneity will persist because a) P), 
equilibrium with the bath liquid at this point, the | the low thermal conductivity of the liquid. Ay hel 
pressure indicated on the manometers will be that | thermometer below the surface will therefore be a nes 
corresponding to the temoerature at the liquid sur- | a lower temperature than that corresponding tf} gy 
face and not at the peramagnetic salt (but see below). | the measured value of p. The paramagnetic salts} yj 

As the temperature of A changed with changing | a much better heat conductor than liquid helium} — gy, 
bath pressure, the variation in susceptibility was | but, even so, a very small heat influx (if all passe) jj), 
detected as a change in the bridge setting, ».  Sinee | through the salt) will suffice to set up a differential con 


the approximate calibration formula of the thermom- | of several millidegrees across a 1-in. sphere. (4 sms 
eter was 7'=34/n and the precision of bridge setting | rough calculation gives: at 4.2° K, 5 x 107* watt pe 
was An=0.001, temperature changes could be | millidegree; at 2.2° K, 107* watt per millidegree.[ + 
observed as follows: At 4° K, 5107 deg; at | The employment of a vacuum jacket in order t 
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avoid a cold spot on tube F (fig. 1), was found to do 
more jiarm than good; heat flow along F from above 
js prevented from entering the bath-liquid, and the 
conduction heat leak into the vapor-pressure bulb 
js greatly intensified. A strong manifestation of 
this elfect was observed in the first apparatus, and it 
was found necessary to admit “exchange gas’’ into 
the erstwhile vacuum space, G, in order to minimize 
this heat leak. [In facet, the “cold-spot hypothesis” 
appears, upon closer examination, to be invalid for 
liquid-helium vapor-pressure thermometry in a glass 
apparatus. This was not appreciated, however, at 
the time of designing the first apparatus. The 
hypothesis supposes that a pressure drop exists 
along the tube corresponding to the full hydro- 
static head, which is, for example, 1 mm of Hg for 
a 10-cm depth. (This is equivalent to 1.5 millide- 
grees at 4° K, 3.7 millidegrees at 3° K.) 
up the tube, which, for the maintenance of the sup- 
posed conditions, must cool to the cold-spot tempera- 
ture and then return downward by convection or 
condense and flow back as liquid. A rough caleula- 
tion shows at once that in a tube of the size used 
(l-em i. d.) the required heat transfer could not 
possibly take place through the wall and the process 
must be self-stifling.| 

The presence of gas oscillations in the tube, F, 
was sought for by connecting a small oil manometer 
of 2-mm bore between F and the bath. None 
could be detected with the arm connected to the 
bath closed off. (With this arm open, the oil 
meniscus oscillated with a period of some 15 see and 
an amplitude of about 1 mm.) The same manom- 
eter could be used as a direct check on the hydro- 
static-head effect. In the first apparatus, the hydro- 
static head as measured in this way was always 
considerably higher than that computed from the 
level of the liquid in the bath. With the second 
apparatus, the two values always checked much 
more closely, which suggests that a heat leak was 
present in the first apparatus. 

The second apparatus was accordingly — built 
without the vacuum shield, G (fig. 1, a) and with 
smaller holes in the radiation shields (which were 
also increased in number and more closely spaced). 
Enough liquid was condensed to reach, initially, 
the equator of the sphere, A, to make better use of 
the thermal conductivity of the salt as opposed 
to that of liquid helium. 

The manometers were of sufficiently large bore to 
obviate the necessity for meniscus-height corrections. 
Thermomolecular pressure effects, and the reflux of 
helium due to film flow in the He mn region, were 
negligible for the l-cm-diameter tube used. Pres- 
sures were measured to 0.01 mm by means of a 
Wild cathetometer and the readings corrected to 
standard gravity and 20° C. Pressures were stabi- 
lized through simultaneous adjustment of a_fine- 
control pumping valve and the current through a 
small heating coil in the bottom of the cryostat. 





ndebted to J, R. Clement for suggesting this procedure. 


Under | 
such a Ap, there would be a very large flow of vapor | 


The bridge reading, n, is a linear function of the 
susceptibility of the paramagnetic salt specimen, 
x, which varies as the inverse of the absolute tem- 
perature (Curie law). Departures from the Curie 
law due to the crystalline field splitting of the ground 
state spin-quadruplet are negligible down to 1.3° 


K, the lower limit of measurement [14]. The 
calibration formula actually has the form 
n—B=A/T, (1) 


where A is proportional to the Curie constant of the 
material, and B is a second constant, equivalent to 
the bridge balance value at infinite temperature. 
In order to examine the over-all consistency of a 
given p-T relation the two constants A and B may 
be determined by plotting n as a function of 1/7 and 
fitting a straight line to the data. Inserting these 
values into eq (1), magnetic temperatures, 7',, are 
then computed for each measured value of n, and 
the differences A7 = T— T,, provide a measure of the 
“over-all validity” of the given p-7 relation. [As 
the latter may be a priori in error in any part of the 
temperature range investigated, the correct values 
of A and # are, within limits, infinitely variable and 
the final choice correspondingly arbitrary. The two 
constants may be fixed by considering the scale to 
be correct at any two chosen points, preferably at 
opposite ends of the temperature interval (ef. Erick- 
son and Roberts [6]).] Values of A7/7,, are then 
plotted against 7, for a more sensitive check of the 
quality of the above fit: for let us suppose that the 
derived A7’s are entirely due to incorrect choices of 


the values for A and B. Then 
es SC oT. 
AT=T-—T,, “54 S44 oR AB 
= T,, +48 . & (2) 
A A 
and Al _AA 4B a (2a) 


i &° # 


Bearing in mind that this type of plot will tend to 
exaggerate deviations at the lowest temperatures, 
the best straight line may now be drawn and the 
originally chosen values of A and B modified accord- 
ingly (see footnote 5). By such a procedure it was 
found possible to assign final values to A and B with 
confidence that any alternative choices that could 
be readily permitted would only affect the values of 
AT in a minor way and would not invalidate the gen- 
eral conclusions as to the over-all correctness of the 
scale. [This procedure becomes the more accept- 
able, the better is the p-7 relation under examination. 
In the case of the 1948 Seale, for example, the errors 
are so large that the choice of B would have to be 
supported by an independent determination of this 
quantity from measurements at high temperatures; 
this was not possible in the present apparatus where 
the measuring coils were immersed in liquid helium.] 
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TABLE 1. 


( orresponding values of pressure in millimeters of mercury, p (correcte d to 20° C 





dial-setting, n (corrected for nonlinearity in decade scale) 


June 7, 1055 June S, 1055 


THs. OO 11. 22 748. 67 11. lei Ss. 04 
L.@7 11. 178 414. 42 O STH os, 12 
27.40 10. SONS 11.000 * SOs si) 
Vis. 12 Ww 4 7. SO4 . Sao O08. 56 
Mal. 1S 9. 4457 23.473 2. 182 AOL. 32 
24 04 8. 7H 16. YS7 1. 24 W235 
! bb | 7.47% 13. 827 0 AS 173. 72 
71.02 7s l sl s SAO 
70.71 7H 8. 52H 7u S7_ SI 
ww 2? 478 SSH TAS 14. 47 
25. Si 2.451 6. 415 l 4 4.774 

“uw 73 

2. 229 4. 927 

1. 26 6. 76 


4. Results 


The experimental data, 1. e., corresponding values 
of pressure in millimeters of mercury (corrected to 
20° C and standard gravity) and bridge reading 
(corrected for nonlinearity), are given in table 1. 


4.1. Comparison With the Clement Equaticn 


Four runs on different days were made with the 
first apparatus and the results could be harmonized, 
using slightly different values for the constants A 
and B, with a quite small scatter. This is shown in 
fig. 2,a.° 

The signal feature of these results is a large posi- 
tive deviation hump (A7’~5 millidegrees) in the 
middle of the Het region. The deviations are vers 
small in the He 1 region (less than 1 millidegree and 
achieve significant) negative values above 4°K. 
From what has been said in the preceding section, 
the good fit in the He im region could be fortuitous, 
and one might, for example, reassess the thermom- 
eter calibration to give a zero deviation in the 
region of the boiling point and again somewhere in 
the middle, sav, of the Hem region. This has been 
done in figure 2, b. The main effect is to exaggerate 
the 3°K “hump” while only slightly affecting the 
points below the lambda point. 

On either assessment, however, the suggestion 
remains that the “Clement temperature” (7) is 
too high by several millidegrees in the neighborhood 
of 3°K. Equally well, of course, the measured vapor 
pressures could be too high here due to errors of 
measurement, such as would be occasioned by @X- 
traneous heat influx to the bulb. <All the data of 
figure 2 were obtained with exchange gas introduced 
into the jacket, G (fig. 1, a). Without exchange gas 
the 3°K hump became several times larger, and this 
finds a ready explanation in terms of heat conducted 
down tube F (fig. 1, a), as discussed in section 3. 
By the same token, the hump of figure 2 might be 
due to an ineradicable heating effect, and that it is 
partly so is supported by the data obtained with 
the second apparatus, shown in figure 3. 


rhe results as presented in fi 2a and 3 were actually obtained using t} 
came 1. F# values as finally chosen for the corresponding Van Dijik-Durieu 
evaluation (fi fand 5, respect vely It »> happens that these valus pro I; 
Ss good neral fit a inv other 
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Only one experiment was carried out with this 
appar: tus, and the data are accordingly somewhat 
spars Che thermometer calibration was quite 
different from the previous case, i. e., one was then 
operating upon a different section of the bridge 
windings, Which is an advantage in checking the in- 


ternal consistency of the bridge. The results are 
verv similar to those of figure 2, b, but the hump in 
the Her region is much less pronounced (see foot- 
note © 


4.2. Comparison With the Van Dijk-Durieux Table 


Recently Van Dijk and Durieux (see footnote 3) 
have recalculated the p-7 relation for helium, making 
a reassessment of the best available thermodynamic 
data. Their results are not in exact accord with the 
Clement equation and, in facet, deviate from it in 
the Het region after the manner of the points in 
figure 3, with a ATi. of 3.6 millidegrees at 2.9°K. 
In the Hea region, A7(=Te,—Typp) increases 
steadily from 0.001 deg at 1.9°K to + 0.0016 
deg al L.O°K, 

In consequence, the magnetic thermometer data 
of this report can be fitted somewhat more closely 
to the Van Dijk-Durieux table than to the Clement 
equation. Figure 3 has been obtained in a rather 
subjective manner, with an eve to the general 
features of a plot of 7e,— 7ypp versus 7. 
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corre sponding lo June data of 
(Compare fig. 3. 
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The latter [not shown] shows zero deviation at 
1.7°K and a minimum deviation at about 3.8°K; 
these features have been used as a guide in obtaining 
the calibration curves from which figure 4 is de- 
rived. It can be seen that a two-point tie-down is 
forced on one here, as it is impossible to choose values 
of A and B&B (eq (1)) that will give effectively zero 
AT over a significantly large temperature interval 
without causing unduly large deviations elsewhere. 
Thus one may keep the magnitude of A7(= Typp 
T small throughout the entire range of measure- 
ment but then one also obtains a systematic deviation 
in the He 1 region. 

The August data are plotted in figure 5, which is to 
be compared with figure 3. Figure 5 was obtained 
as the best fit throughout the entire range of measure- 
ment, and it, too, has the feature of practically zero 
deviation at 1.7° and 3.8° K. The magnitude of 
AT now lies within +1.5 millidegrees throughout the 
range of measurement. 


5. Discussion 


The two sets of data, viz., the June results with 
the first apparatus and the August results with the 
second, do not agree quantitatively but do show 
qualitatively the same behavior. It is felt that the 
August data are the more valid (smaller heat leaks 
into vapor-pressure bulb), and these agree very 
closely with the Van Dijk-Durieux calculated p-7 
table. Both Van Dijk and Clement have recently 
made significant modifications to their respective 
p-T tabulations. Ina subsequent report the present 
authors will reevaluate the above data together with 
the results of a new series of measurements, 

None of the present measurements can be brought 
into very close accord with the Clement equation, 
in contrast to the results of Erickson and Roberts 
\6|. The greater part of the latter data was obtained 
from bath-pressures plus hydrostatic-head correc- 
tion, which might account for the discrepancy, since 
the conditions necessary for the exact validity of 
applying such a correction are impossible to achieve 
in practice. 

The significant differences among the data of 
Erickson and Roberts and the first and second series 
reported here point up a major problem in vapor- 
pressure thermometry in the liquid-helium region, 
viz., the difficulty of reproducing results from one 
apparatus to another and of checking any p-7' rela- 
tion to the desirable accuracy of 1 millidegree. From 
the practical point of view, one desires to obtain 
accurate values of absolute temperature from simply 
performed measurements of vapor pressure, the 
latter being an essentially minor part of any given 
investigation. It is not entirely unreasonable to 
seek a solution in the direction of making rigid 
stipulations concerning the technique of vapor- 
pressure measurement and to evolve a “practical 
p-T relation” that is different from that calculated 
on a thermodynamic basis. 

On the other hand, it is true that the differences 
between the Clement equation and the Van Dijk- 
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Durieux table are never greater than a few milli- 
Allowing a 2-millidegree precision for the 
latter table (as its authors claim), a “Clement tem- 
perature” should have a maximum deviation from 
the thermodynamic temperature of 5.6 millidegrees 
at 2.9° K, a probable error that is tolerable for many 
of the common research problems in this region. 
demanding a knowledge of d7/dp 


degrees 


Experiments 


e. g., specific-heat determinations) are, of course, 
most sensitive to deficiencies in’ the temperature 
scale. 

The authors are greatly indebted to J.) R. 


Clement of the U.S. Naval Research Laboratory for 
many stimulating discussions and helpful sugges- 
tions. It is a pleasure to acknowledge the further 
benefit derived from general discussions of the prob- 
lem with R. A. Erickson (Ohio State University), 
M. P. Garfunkel (Westinghouse), P. H. Keesom 
(Purdue University), W. E. Keller (Los Alamos), 
lL. D. Roberts (Oak Ridge), C. A. Swenson (lowa 
State College); and with H. van Dijk, D. de Klerk, 
and M. P. Durieux (Kamerlingh Onnes Laboratory). 
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Calculation of Thermodynamic Functions for 
Polyatomic Molecules 


Harold W. Woolley 


Formulas are given from which the thermodynamic functions may be obtained for 
polyatomic molecules not exhibiting special phenomena, such as internal rotation. The 
effects of first and second order anharmonicities and rotational-vibrational interaction are 
treated in detail, as also are the angular momentum effects of doubly degenerate vibrations 
in linear molecules. This paper presents formulas for the corrections carried to higher order 
than has hitherto been available. 


l. Introduction 


There is occasionally need for accurate thermodynamic functions for polyatomic molecules. 
A definite calculating procedure exists for approximate calculations via the rigid-rotator har- 
monic oscillator treatment, with small corrections for effects of anharmonicity, stretching, etc. 

This paper presents formulas for the corrections carried to higher order than has hitherto 
been available. Calculations can thus be made via the approximation formulas with accuracy 
probably comparable with that of the direct summation method using extrapolated formulas 
for the energy levels, aside from effects related to vibrational and rotational cut-off, which have 
been omitted. 


2. Nonlinear Polyatomic Molecules 


The partition function for the internal molecular energy states may be written as 


= OrOren .o. Qe. (1) 


For nonlinear molecules, 
Qre= (1 y) y (/ Ap Polo) (k They (2) 


is the rigid rotator partition function for the lowest vibrational energy state, where Ao, 2, and 
(, are the rotational constants for the ground vibrational state, and y is the symmetry number. 
Quo. is the partition function for vibrational levels of harmonic oscillators. Qs, is a correction 
factor including (1) a factor 1+ p7. . . for rotational stretching, as shown by Wilson [1]? and 
(2) a factor 1+ 6, 7+ 6,7? for a low temperature quantum correction for rotation as shown by 
Stripp and Kirkwood [2].° Q, is a factor including rotational-vibrational interaction and the 
anharmonicity corrections, now to be treated in detail. 

The rotational-vibrational interaction in its simple form arises from the vibrational de- 
pendence of the analogue of eq (2) for higher vibrational states. Before summation over 
Vibrational quantum numbers, the dependence is introduced by including the factor 

P, exp|La v,+3a,07+34,00;+- 354a,,0°+ Sa, 070,454, ,00 ry), (33) 
equivalent to a form introduced by Benedict [3].!. With his formulation, the constants in the 
factor A, come from sums of corresponding constants for .1,, ?,, and C,. A summation pro- 
cedure adequate to treat this with the anharmonicity effects follows, extending the latter 
bevond the linear terms of Stockmaver, Kavanagh, and Mickley [4]. 


Related to a thesis supmitted to the Graduate School of the University of Michigan in partial fulfillment of the requirements for the Ph. D. 
le gree in Physies 
Figures in brackets indicate the literature references at the end of this paper. 
Expressions representing the 7? rotational stretching effect and the 7! low temperature correction as well as a low temperature correction 
rsymmetry effects for different classes of rotational levels for the symmetric rotator are included in the thesis 
‘An equivalent representation in slightly different form is given in the thesis 
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For a given vibrational level, specified by vibrational quantum numbers ¢,, 1 <i <n, the 
statistical weight is the product of the n quantities 


p vr, +d,—1)!/[ef(d,;—1)!], (4 


where d, is the degeneracy associated with the 7th fundamental frequency. 
The summation over the vibrational levels may be indicated as 


t 


i) OrQrerR exp| he, kT \Wp (5 
ia i 


with the vibrational energy in wave numbers given by 
G, = 0.77 oF PA 1) SJ rr ZY PAG, L)(r 2) ; ay rie le. Ly ue Te (()) 


The relations between these constants and the more usual ones are indicated in appendix 1. 
The Boltzmann factor is separated into two factors, one for harmonic oscillators fitting 
the lowest levels and a second giving the alteration due to anharmonicity 


exp[—heG,/kT]=exp|—heSo,r,/k Tlexp[—he(G,— wr) /kT7). (7 


(J... the product of the second factor on the right hand side of (7) with P,, is expanded into 


a power series Whose terms are polynomials in the ¢’s. Defining a function 


{(Z)=r,!/(v,—Z)! (S) 
it can be shown that 
Rah han (Ss 9) 
where u,~—hew, kT, and that 
Dut (Z) pit “."% —[(d,+-Z—1)!'/(d 1 )!]-[e-4 l—e lied | (10) 


The factor Q,, can be expressed as a sum of terms, each linear in whatever f,'s it contains.’ 
The results of multiple summation over all ¢,’s follow by application of eq (10). If the right- 
hand side of eq (10) be designated by f*(Z), one has 2,.Q,Q.. = Qu o.Q-. Where Q, = exp [— 24,7, 
Quo. —=,Q, and where Q, is the same function of f/*(Z) as Y,. was of f,(Z). Thus eq (1) 
has been evaluated, Qy 6, being given by 


(0.0 I] l f q; (11) 


There is some advantage in using the logarithm of the partition function from this point 
on because of analytic cancellations of some higher order terms. The terms removed are cross 
products without linkage in the subseripts, sueh as sjjsy:, £)%),, and ay)”, for example. 

The contribution to In Q, due tov, and v,, terms, with r,=« and s l—< Is 


(he kT)>° rd f(d,+1)rs* r,ddrarss ; the 2 } 2x: dj(d,- L)r7tl +-P2i(d, l)r |x 








—~Pddrrilady dor |s?s? +277. ddyr rr s?s 8,2 42,,0,(d + 1)d rrjs's (12) 
with terms of higher order for this and other contributions listed in appendix 1. At high 
temperature the first term varies as 7 and the second as 7°. The terms in the appendix are 


higher than 7° in principal dependence at high temperature. 


See appendix 2 
Relations giving the simpler polynomials in tert the f’s will be found in appendix 
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Similarly, the contribution of higher order anharmonic terms is 


(he/kT) S30 f ye (d+ 1) (d+ 2 )r3s83 + yd (d+) d jrPrjs?s + y jd ddr ir rs 88x}. (13) 


(he contribution due to rotational-vibrational interaction is 


Do iadyrs;+a,dr (itd ys? +a, dy ryss,4- abd rs? (14) 


as far as a’s alone are concerned, while their mixing with anharmonicities gives 
(he/kT)>) { 2a yr, (dy 4-1 )r?83 + a ry dr rys?s;} (15) 


In all of these (multiple) summations, if the terms are completely symmetric in the indices 
that are present, such as 7, 7, ete., the symmetric indices are to progress in size as with << j<k, 
and so forth. For terms not completely symmetric in the indices, the indices not symmetric 
are simply not to be equal, thus 1#j#k, and so forth.’ 


3. Linear Molecules 


There is an additional kind of vibrational-rotational interaction, concerning vibrational 
angular momentum, which is important for linear molecules. When doubly degenerate vibra- 


tions occur, the vibrational energy, eq (6), contains an additional term g=4q,J/d;. For 
&. | * J 17 J 
ie 


linear molecules, in addition, there is dependence of the rotational energy upon / where /= X/,, 


viving for the two together 


F(J, l)=9+ B,[J(J+1)—P]—D,|J (J+ 1)— PY, (16) 





where J=/1,/+1,/+2 .. ., with double degeneracy for 1 >0. 
The sum over the rotational levels for J>‘/, follows from a half integer summation formula 


of the Euler-Maclaurin type [5] 


>> f(n) | f(n)dn + (1/24) f" (ny — 3) — (7/5760) 7" (ny — 3) +- (31/967680)f © (ny —3) 2. ., (17) 


viving 





>) (24+ 1)exp heBlLJ(J+-1)—PVkT=k TT /heB+1/38+ ... (18) 
J=] 
and 
Q= 3 (2S + De MF T= (1/ yk T /heB + 1/3 +-2(k T/heB)2D/B) ee", (19) 
J 


The exponential in g is expanded as a power series and the summations made with 7, increasing 
by steps of two from /, to ©, then summing /; from — o to +o. If gq is extended to include 
possible contributions to In Q varving as 7° at high temperature, while retaining the g,; term 


needed for nonlinear molecules, as in 


a= Dog? + D9 d+ Dog udPeit Do iake:, (20) 


t 


the contribution of q to In . becomes 


2(he/k TYE | gyri? geri +39 sa gid et S28_— 29,0 7288 } + (he/kT Pdf g?; (r+ 87? 


3\ 4 972 29 >. on 2» ‘ —s » \e3o. ) 91) 
+r? st + 2g? rr ysis? 4g tr? (4+ 3r)st+ 29 rar re +rpsis, }. (21) 
For a case with 3 frequencies, for example, which is not completely symmetric in the indices would include r),27)3 and also rggr But 
; has complete symmetry in the indices, and rysri,7), is sufficient without rygry.7y) also 
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Constants of the latter two types in eq (20) are not experimentally available at present. For » 
linear molecule, when there is only one doubly degenerate frequency, with /=/,, the — P,/ 
term of eq (16) is commonly absorbed into the goo/3 term. In that case, g;, of eq (21) repre 
sents gx +P, of the spectroscopist. 

The free energy function for a mole of a substance in the standard ideal gas state is given by 


(F°— fy)/RT=In Q+ (5/2) In T+ (3/2)In M-+In (22/N)°?K5?/h3 Po, 
P, being atmospheric pressure, \/ the molecular weight, and the others known physical con- 
stants. The other functions follow by differentiation as given in appendix 5. 


A tabulation by Pennington and Kobe [6] may also be used for convenient evaluation of 
the more important of the contribution here treated. 
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5. Appendices 
5.1. Appendix 1 


The relation of the constants of eq (6) with those in 


G SJa,(r,+d 2)+ S57 (vy +d 2)°+ QF (v,+d,/2)(r,+-d 2)4+ SO iiileit+d,/2) 
aay / d,/2)2(v,+d 2)+ Day .(v,+d,/2)(v,+d,/2)(r,+-d,/2) 

is 
w=0,+ (d,+1)7,,4 > 3,2 i+ (1/4)(3d? + 6d,- 4)7 

| } DO )(d t Ld jy ; D5 (dP /A) Hs; , 22 25 (did, 4)7 - 

I Laat (3) (3d +O), + So(d 2)y 
r Tiit(d l)y t-(d 4 1)y + old 2) 
7] Vii Y icy ANA Ys je — Yue: 


5.2. Appendix 2 


To obtain eq (9) and (10), begin with the identity 


SS 27°=(1—2z7"') 


—_—— 
’ “ 


Multiply by #77! and differentiate » times with respect to x, getting 


. , ! 
al 1)’ (vw-y Z)' , s—148-—e d’ 
r=0 \r 7)! dr” 
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For the second member of the identity, let s—1=t, giving 


d" [(¢ it 1)7¢ | or d” <* : Z| {2 P=i 
dt” dt" | <7 (Z—r)!r! 


, ‘ : ; ig : P Sa d” ; 
For »>Z—1, the only contribution is from the one term with r=Z, giving yal or (—1 )"n'! 
( 
r Z+n 
s—1)-'"". Multiplying both members of the identity by (—1)” » one has 


(n—Z)! 


, (vt+n—Z)la-" n! 


S YY 
<j (v—Z)!(n—Z)! (n—Z)! 


*(1 


Inserting e~“: for s+, one may obtain eq (10) by using n=Z-+d,—1, with eq (9) as the special 
case given by Z=0. 


5.3. Appendix 3 


Functions of ¢ in terms of (f(s) =e!l/(r—s)!: 


e=f(1) 
m—f(2)4+fQ1) 
m—=f(3)4+-3f(2) +f) 
v'=f(4)+6f(3) 4+ 7f(2)+f0) 
v= f(5)+-10f(4) +-25f(3) + 15f(2) +f) 


v= (6) + 15f(5) + 65f(4) + 90f(3) + 31f(2) +01) 


we—r=—f(3)+2f(2) 


vi—v—f(4)+ 5f(3)+-4f(2) 


*— 2 +7? = (4) +4/(3) +2f(2) 
py —2r' +r =—f(5) 4+ 8f(4) + 14f(3)+-4/(2). 


5.4. Appendix 4 
Additional terms in the formulas given earlier are as follows: 
Equation (12): 
(1/6) (he/k TY 33 { 443,d(d 4-1)? [1 + (Sd, + 12)r, 4 (Sa? 4+- 22d 4-15)? + (20? +-4d 4-2)? ]s? 


Had dri (8d, 4+ 1)r,4- (8d)4+ 1) +r? +r + (6d dj 4+-3d,4-3d;+ Iraj +@rr,; 


Pr r?|s3s3 + 2a? rd (d,4-1)dyPr {1+ (3d,+-4)r,4+-(d,+1 \r?|s?s ' 6a, a7,d (d, t 1 )dr?r [2 
(2d,4+-1)r,4-2d,r,4+-dyryr|sts? + 3.7 ad djdyr ry + (2d, + Ir, 4-d jr +d rj|sis?s, 


+ 249 prj (dy + 1d (d+ 1) r2r3s3s3 + 12a yr yt ed (d+ 1) dd yr2r jr, (2 +7 )stsys; 


j 


12r 2,2 d(d,+1 \d dyrer rpSisesy + On wad dddiryr irri t+ 7 SPS 8,8, 


i 


Oa yt yt pd dydydur rr i922 sp8. + 6a yt yt pd d dy yr yr +d rt dr + dyn, 828785 
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1) (d,4-2)d rer js3s; 
Sayer dd ydyr ir pr Pi8 S881 
1)d?r,{1 


Equation (13): 
(he kT)>3 Sti id (d, 1 1)(d, T 2)(d, +-3)rts* ry (d, t 
t 2; d(d, ; 1 )d,(d, : L )reresis? | 


>, A (d | Vdd pr 7 iP p82 8 58) 
(1/2)(he/kT)S3 { 122, yud (dit 1) (di 4-2) [1 + (di 4+ 1) rs? +42, 45d (d 4 
1 )d (d+ 1 )r?r2s8s? +40 yind (d+ 1)d dyr2rr,s3s)s; 
dr, + 2d jr,|s3s? 


22; Yi (d; T 
22 Yurd daddy x ris? 


22 Yer djd,(d, 


I)dy2r,l24 


Lr rss; 


4r,yj (dA 


2(d, T l Ir |sts, 
t+ ba Vid (d, t 
42,Yind(d,4 


27 Yind ddr wry (A 


1)(d,+2)d rer sts, 
<? ? 
A 
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